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ABSTRACT 


An  analysis  and  computational  procedure  has  been  developed 
tor  predicting  the  natural  vibration  modes  and  frequencies  of  rotor  blade 
and  hub  systems.  This  method  is  capable  of  treating  rotor  systems 
which  are  01  a  general  configuration  and  of  including  configuration 
variables  which  are  normally  neglected. 

An  experimental  apparatus  and  a  testing  technique  have  been 
developed  for  obtaining  the  experimental  data  necessary  to  verify  the 
analytical  procedure.  This  apparatus  and  the  technique  are  capable  of 
determining  the  natural  vibration  modes  and  frequencies  of  two-bladed 
rotor  systems  of  u  very  general  configuration  while  operating  in  the 
absence  of  aerodynamic  loads.  An  initial  test  program  with  the  appara¬ 
tus  was  conducted  to  determine  the  adequacy  of  the  equipment.  The 
tests  not  only  demonstrated  the  capability  of  the  equipment  and  test 
procedure  but  also  revealed  a  problem  with  extrinsic  excitation  of  the 
test  rotor  at  frequencies  other  than  the  shaking  frequency.  Recom¬ 
mendations  are  made  for  modifications  to  the  system  and  the  techniques 
in  order  to  reduce  the  levels  of  these  extraneous  excitations  and  to 
minimize  their  effects. 

The  computed  and  measured  spanwise  moment  distributions 
and  the  corresponding  natural  frequencies  are  presented  for  six  of  the 
natural  vibration  modes  of  the  test  rotor  at  several  rotational  speeds. 
The  theoretical  and  experimental  results  are  compared,  and  their 
differences  are  discussed. 

The  model  rotor  used  in  these  initial  tests  was  intended  to  be 
so  simple  as  to  be  representable  by  the  classical  uniform  uncoupled 
beam.  Because  the  theory  for  this  classical  configuration  is  well 
established,  the  differences  between  the  measured  and  computed  results 
were  to  be  assigned  to  the  experiment.  However,  the  results  of  the 
tests  indicate  that  the  test  rotor  could  not  be  adequately  represented  by 
the  classical  uniform  uncoupled  beam,  especially  in  the  higher  frequency 
natural  vibration  modes. 
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INTRODUCTION 


The  structural  dynamic  characteristics  of  helicopter  rotors 
have  a  significant  influence  on  the  vibratory  loadings  experienced  by  the 
rotor,  the  control  system,  the  fuselage,  and  fuselage  contents.  Thus, the 
ability  of  the  designer  to  predict  adequately  these  dynamic  characteristics 
is  essential  for  minimization  of  the  vibratory  loads  experienced  by  the 
aircraft  and  its  components. 

The  many  configuration  variables  and  physical  parameters  of 
the  blade  and  rotor  system  apparently  provide  the  designer  with  a  great 
latitude  to  the  layout  of  a  new  rotor  for  satisfying  the  various  perform¬ 
ance  and  structural  integrity  requirements.  In  practice,  however,  the 
uncertainties  involved  in  the  predictions  of  the  aerodynamic  loads  and 
structural  characteristics  encourage  an  extremely  conservative  approach 
in  which  maximum  use  is  made  of  features  which  have  previously  proved 
to  be  satisfactory  to  the  particular  manufacturer. 

Extremely  conservative  practices  with  regard  to  rotor  design 
have  been  justified  on  the  grounds  that  changes  to  the  geometry  or  elastic 
characteristics  of  the  rotor  would  change  the  aerodynamic  forces  in  an 
unpredictable  way.  Consequently,  radical  departures  from  established 
successful  rotor  systems  have  been  avoided.  A  fair  amount  of  evidence 
suggests  that  this  judgment  has  been  correct.  Recent  developments  in 
the  aerodynamic  theory  of  rotor  blades  and  the  beginning  of  an  under¬ 
standing  of  tip  effects  suggest  that  aerodynamic  forces  are  now 
sufficiently  predictable  to  encourage  a  more  precise  approach  to  blade 
layout.  There  is  a  corresponding  need,  however,  for  the  establishment 
of  a  validated  structural  dynamics  theory  that  will  permit  the  combining 
of  the  new  developments  in  aerodynamics  with  reliable  structural  and 
stability  analyses  to  arrive  at  rotor  designs  that  satisfy  a  number  of 
simultaneously  imposed  constraints.  The  effort  reported  herein  is  the 
first  phase  of  a  program  initiated  to  develop  and  validate  a  reliable 
method  for  predicting  the  structural  dynamics  of  rotor  systems. 


At  the  present  time,  blade  articulation  is  largely  a  matter  of 
individual  company  practice,  and  the  particular  methods  of  attachment 
and  hinge  arrangement  have  evolved  from  cut-  and-try  procedures.  The 
distribution  of  structural  material  reflects  stress-level  estimates,  and 
the  resulting  blade  designs  tend  to  have  some  variations  in  the  mass  and 
elastic  characteristics  which  are  generally  ignored  in  current  vibration 
analyses.  For  example,  the  use  of  bonded  doubler  skins,  particularly  in 
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location.  Also,  point  forces,  such  as  those  introduced  through  a  pitch- 
control  horn,  usually  are  neglected  in  the  estimation  of  rotor  and  blade 
dynamic  characteristics.  Finally,  there  are  virtually  no  quantitative 
data  which  would  permit  the  evaluation  of  structural  dynamic  analysis 
methods  of  rotor  and  blades  under  the  ideal  condition  assumed  in  their 
development- -that  is,  in  the  absence  of  all  aerodynamic  loads. 
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One  purpose  ol‘  this  effort,  then,  was  to  develop  an  analytical 
procedure  for  predicting  the  natural  vibration  modes  and  frequencies  of 
representative  helicopter  rotor  systems  (operating  in  the  absence  of 
aerodynamic  loads)  ot  general  configuration.  The  other  purpose  of  this 
effort  was  to  develop  an  experimental  apparatus  and  a  test  techniq  e  for 
obtaining  data  to  verify  the  analytical  method. 

The  analytical  procedure  which  has  been  developed  is  a 
generalization  of  the  associated  matrix  analysis  technique  (reported, 
for  example,  in  References  1  and  2).  The  theoretical  effort  was  largely 
devoted  to  the  derivation  of  expressions  for  the  inertia  loadings  and 
elastic  restoring  moments  (used  to  define  the  elements  of  the  mass  and 
elastic  matrices)  for  a  rotor  and  blade  of  a  very  general  configuration. 
The  analytical  procedure  includes  the  effects  of  the  following  rotor 
parameters : 

1.  Nonunilorm  mass  and  torsional  inertia  distributions. 

2.  Nonuniform  flapwise,  edgewise,  and  torsional  stiffness 
distributions. 

3.  Built-in  twist. 

4.  Mean  bent  shape. 

5.  Noncoincident  and  nonstraight  cross-section  c,  g. , 
centroid,  and  elastic  axes  which  do  not  pass  through 
the  rotor  axis  of  rotation. 

6.  Details  of  articulation  including: 

a.  Root  fixities. 

b.  Radial  and  chordwise  hinge  offsets. 

c.  Hinge  inclinations. 

a.  Pitch-horn  location  (point  forces), 
e.  Pitch-axis  offset  and  inclination. 


The  experimental  effort  was  largely  devoted  to  the  design, 
fabrication,  and  installation  in  a  vacuum  tank  of  a  rotating  shaker 
system  and  rotor  drive  and  support  system.  The  initial  experimental 
investigation  reported  herein  was  intended  primarily  to  establish  the 
adequacy  of  the  apparatus.  This  was  to  be  accomplished  by  measuring 
the  natural  vibration  modes  and  frequencies  of  a  model  rotor  system 
intended  to  be  simple  enough  to  be  adequately  represented  by  the  classi¬ 
cal  uniform  uncoupled  beam.  Because  the  theory  for  the  dynamic 
response  of  this  classical  configuration  is  sufficiently  well  established, 
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any  differences  between  the  measured  and  computed  results  could  pre¬ 
sumably  be  assigned  to  the  experiment.  The  results  of  these  tests  and 
comparisons  with  the  simple  theory  are  presented  and  discussed. 
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THEORETICAL  DEVELOPMENT 


INTRODUCTION 

It  is  the  purpose  of  this  analysis  to  develop  a  method  for  pre¬ 
dicting  the  natural  vibration  modes  and  frequencies  (in  the  absence  of 
aerodynamic  loads)  of  rotor  blades  which  are  of  a  general  configuration. 
The  analysis  will  allow  for  the  combined  elastic  and  rigid  body  motions 
of  the  blades,  arbitrary  offsets  of  the  root  hinges,  and  noncoincidence 
of  the  blade  twist  axis,  elastic  axis,  and  centroidal  axis.  Furthermore, 
these  three  axes  need  not  be  straight  and  need  not  pass  through  the 
center  of  rotation.  The  elastic  description  of  the  blade  includes  flapwise 
bending,  chordwise  bending  and  blade  torsion,  while  the  rigid  body 
degrees  of  freedom  are  flapping,  lead-lag,  and  rigid  body  pitching. 

The  basic  nomenclature  used  in  the  analysis  are  defined  in  the 
list  of  symbols. 

The  second  section  is  concerned  with  the  elastic  behavior  of  a 
twisted  blade.and  formulas  are  developed  for  the  elastic  restoring 
moments  as  functions  of  blade  deformations.  Expressions  are  also  given 
for  the  deflection  of  a  mass  point  in  the  blade  relative  to  a  blade  refer¬ 
ence  axis  system  which  rotates  with  the  blade.  The  treatment  of  the 
elastic  behavior  of  the  blade  has  been  based,  to  a  large  extent,  on  con¬ 
cepts  developed  by  Houbolt  and  Brooks  (Reference  3)  in  considering  a 
simpler  blade  representation  without  rigid  body  degrees  of  freedom.  In 
particular,  the  axial  stresses  due  to  centrifugal  forces  are  treated  as  if 
the  blade  acted  as  a  collection  of  longitudinal  fibers. 

Expressions  for  the  inertia  forces  acting  on  the  blade  are 
required  in  the  procedure  developed  for  finding  the  blade  dynamic 
motions  and  natural  frequencies.  These  are  obtained  by  finding  the 
acceleration  of  a  differential  mass  in  the  blade,  multiplying  by  the  mass 
to  determine  the  inertia  force,  and  then  carrying  out  the  appropriate 
integrations. 

The  problem  of  determining  inertia  forces  is  complicated 
because  of  the  generality  of  the  configuration  which  has  been  allowed  and 
the  need  for  determining  absolute  accelerations  relative  to  an  inertial 
frame.  The  derivation  of  the  blade  inertial  loadings  is  somewhat  lengthy 
but  is  carried  out  by  a  series  of  straightforward  steps.  Successive  sec¬ 
tions  treat:  (1)  the  motion  of  coordinate  systems  by  a  series  of 
transformations,  (2)  the  derivation  of  the  absolute  acceleration  of  a  mass 
point  of  the  blade  from  the  accelerations  relative  to  moving  coordinates, 
and  (3)  the  integration  of  the  blade  inertial  forces  to  obtain  blade  inertial 
loadings  per  unit  span  referred  to  the  blade  reference  axis  system. 
Finally,  the  inertial  loading  is  transformed  to  moments  about  the  local 
elastic  axis  and  the  major  and  minor  axes  of  the  blade  section.  Equa¬ 
tions  are  then  derived  for  the  total  force  and  moment  components  acting 
on  a  section  of  the  blade. 
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The  analysis  has  been  simplified  by  dropping  terms  in  the 
course  of  the  development  which  would  lead  to  nonlinear  terms  in  the 
oscillatory  variables  in  the  final  equations.  It  has  also  been  assumed 
that  the  mean  values  of  the  nondimensional  structural  deformations  and 
the  mean  flapping,  lead-lag  and  pitching  angles  are  all  small  quantities, 
so  that  their  squares  and  products  can  be  neglected  when  they  occur  in 
the  coefficients  of  the  oscillatory  variables. 

General  differential  equations  for  studying  the  vibrations  of  a 
nonuniform,  twisted,  rotating  helicopter  blade  including  elastic  and  rigid 
body  motions  could  be  obtained  by  combining  the  expressions  for  the  blade 
elastic  behavior  and  the  expressions  for  the  blade  inertial  loadings  and 
acting  bending  and  torsional  moments.  Such  a  development  would  be 
analogous  to  the  one  used  in  Reference  3  .  However,  the  resulting 
differential  equations  would  be  so  complicated  (even  though  linearized) 
that  some  type  of  approximate  method  would  be  required  to  solve  for 
the  vibration  characteristics  of  the  blade. 

A  different  type  of  procedure  was  developed  in  the  present 
investigation  for  determining  natural  frequencies  which  is  an  extension 
of  Targoff's  associated  matrix  technique  (References  1  and  2  ).  In  this 
procedure  the  blade  is  represented  by  a  number  of  spanwise  segments  or 
bays,  and  the  inertial  loadings  on  each  bay  are  assumed  to  be  concen¬ 
trated  and  to  act  at  the  center  of  the  bay.  Each  bay  is  then  divided  into 
two  parts--one  outboard  and  one  inboard  of  the  concentrated  mass.  The 
half  bays  are  treated  as  if  weightless, with  the  concentrated  mass  being 
located  at  the  junction  between  them.  The  elastic  properties  are  assumed 
to  be  constant  within  the  outboard  and  inboard  halves  of  the  bay  (but  not 
necessarily  the  same  on  each).  The  built-in  blade  twist  is  incorporated 
in  the  model  by  permitting  angle  changes  at  the  junction  between  bays. 

The  variables  in  the  vibration  problem  include  the  forces  and 
moments  acting  at  the  ends  of  each  bay;  the  elastic  deformations  at  the 
ends  of  the  bays;  and  the  flapping,  lead-lag,  and  pitching  angles.  The 
variations  of  these  quantities  along  the  span  are  related  by  matrices. 
Expressions  for  the  elements  of  these  matrices  are  derived  in  successive 
sections  of  this  chapter  by  application  of  the  results  for  the  elastic 
behavior  of  the  blade  and  the  blade  loadings  obtained  in  preceding 
sections.  The  transformation  matrix  gives  the  changes  in  the  variables 
at  a  section  where  there  is  an  abrupt  change  in  built-in  twist.  Elastic 
matrices  for  the  inner  and  outer  parts  of  a  bay  give  the  changes  in  the 
variables  across  the  inner  and  outer  sections  of  the  bays  respectively. 
Finally,  a  section  is  devoted  to  the  mass  matrix  which  gives  the  changes 
in  the  variables  from  a  position  just  outboard  to  a  position  just  inboard  of 
a  concentrated  weight  at  the  middle  of  a  bay. 

Relationships  can  be  obtained  between  the  variables  at  the  tip 
and  at  the  root  of  the  blade  by  multiplying  the  matrices  for  all  the  blade 
segments  together.  In  the  final  section  of  the  chapter  a  discussion  is 
given  of  how  a  solution  for  the  vibration  characteristics  of  the  blade  can 
be  obtained  from  these  relationships  and  the  application  of  approoriate 
boundary  conditions. 
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ELASTIC  BEHAVIOR  OF  A  TWISTED  BLADE 


In  t hi ^  > e<  t'.on.  the  elastic  restoring  moment  components  of  a 
v  ms-,  section  ut  a  twisted  blade  are  expressed  in  terms  of  deformations. 

I  he  development  is  made  tollowing  the  principles  of  engineering  beam 
theorv  -innlar  to  the  one  given  bv  Houbolt  and  Brooks  in  Reference  3  . 

I  he  blade  ettevtive  «.  ross  section  may  change  gradually  along  the  span 
b'.t  is  assumed  svmmetru  about  the  chord.  The  axis  of  built-in  tw’ist 
\j  *  ).  the  elastic  venter  (c.r  )  and  the  centroid  are  all  on  the  chord.  In 
the  treatment  ot  the  elastic  behavior,  the  present  analysis  differs  from 
the  analvsis  ot  Houbolt  and  Brooks  principally  in  that  the  axis  of  built-in 
twist  and  the  elastu  axis  ot  the  undeformed  blade  are  not  assumed  to  be 
coincident  arii-1  aliened  with  the  z  ax»s.  Due  to  the  built-in  twist,  the  blade 
angle  d  varies  with  z  and  the  distances  ey,  e^.  ef  and  eA  shown  on 
Figure  1  nwv  also  change  grudunllv  along  the  span. 
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Figure  1.  DIAGRAM  DEFINING  BLADE  STRUCTURAL  AXES 
AND  POINTS  IN  BLADE  SECTION. 


The  f  j  }  reterence  axes  for  the  blade  are  fixed  relative  to  the 
uiacn  root.  Equal. uns  are  first  obtained  for  the  displacement  of  a  general 
point  .n  the  "lade  relative  to  these  axes  assuming  that  plane  sections 
remain  jlai.o  us  the  nladc  deflects.  The  longitudinal  fiber  strains  and 
v  or  responding  stresses  at  any  point  in  a  cross  section  can  then  be  derived 
iron  i no  equations  lor  the  displacements,  and  finally,  expressions  can  be 
o'lt.i.i.rc:  tor  tin-  contribution  ot  these  stresses  to  the  internal  elastic 


Before  the  occurrence  of  any  deformation,  let  there  be  a  point 
in  a  blade  section  at  &  -  x0  with  the  coordinates  y  =  and  y  -  ya  . 
Then,  from  the  sketch  on  Figure  1  , 

y0  =  cy  *(et.7j)cos/S  -  5  S/V  (1) 


f0  =  +(et+7j)&wjB  *  £  cos/3 


When  the  blade  is  under  loading,  let  <£  be  the  torsional  dis¬ 
placement  of  the  section  originally  at  z0 ,  let  u  be  the  displacement  in  the 
x  -  direction  of  the  centroid  of  this  section  and  let  V  and  ur  be  the  dis¬ 
placements  due  to  bending  of  the  elastic  center  of  the  same  section  in  the 
y-  and  ^-directions,  respectively.  Then,  the  original  point  ( x0  ,  y0  ,  ) 

will  have  the  following  new  coordinates : 


x  s  x0  +  «  -  v ' 


j(ef-  eA  +  *{)  ~  C  s,‘n(/3*  <f>)  ■ 


-  ur'  -j  (e^  -  eA  +  rj  )  S]n(j3+  <t>)  +  C  cos(/3  *  <p)  j- 


y  =  ey 


f*- .  A  x  d<p  .  \ 

-J  <er^  *17** °  +  ^ 


+  )  cos(/+  $)-  Z  s m(/$  +  <f>) 


er  +r°  ITT  ^ 

(0  ^0 

*■  {tf  +  lj  )  $m(j3+  <p)  f  £  cos{'/+  <f>) 

where  e^or  e y  is  the  value  of  e y  or  ey  ,  respectively,  of  a  section  origi¬ 
nally  at  .  Here  a  prime  is  used  to  indicate  differentiation  with  respe 
to  x0 ,  such  that 


._/  *  tr  _ »  cC 1ur 

It  -  — - —  ,  uT  -  — - -  >  etc, 

*xe  ct  X  q 
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The  quantities  Ay 


and  A , 


are  defined  as 


Ay  =  p  (<=,  -  ty  >  **• 


and 


and  give  the  linear  displacements  of  the  elastic  axis  at  %  =  xa  arising 
from  torsional  deflection  at  more  inboard  stations.  When  torsional  dis¬ 
placements  are  small,  it  can  be  assumed  that  cos  <f>  =  1  ,  =  and 

v’<p  =  b  -  0  ,  and  the  preceding  expressions  for  z,  y  ,  and  y  become 


**  **  *•'(%  cos/S) 


■  -  ey~eA 


(3) 


y  :y»  +  r  Ay 


(4) 


y >o  *  w  *(^~ey)^  *  Ay 


(5) 


Since,  with  <f>  -  0  at  x0  =  0  , 

AV  e -f/‘  “y*'  <** o 

and 

A  y  =  ey  $  ~  J 'r°  Sy  f  '  ci  XQ  , 

it  is  clear  that 

Ay  =  e y  <j>  and  Ay  -  ft  . 
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Thus, 


%'  =  /  *  u  '  -  V- "  ( yo  -  e  y  -  eA  cos /3 ) 
t  -ey  -  eA  sm/S) 

(6) 

“  ■  ey  “  cos/  *  &A  ft' 5 'n/3) 

“W‘'ty-e}-eA  Sm/*-eA  f'  C05/S) 

S's  So  *  -(fo-e?)?'  -  t'o  t 

(7) 

/*  -ey)  *'  *  (8) 

where  </„'  and  ^  are  evaluated  for  constant  >7  and  £  .  By  differentiating 
(1)  and  (2),  it  is  found  that 

Sj  =  ey  '(?o  ~  e?  )  +  e  't  cos/6 

to  =  e\  +/*'(So-ey)  +  e's  *"7* 
and,  hence,  Equation  (6)  becomes 
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Now, 

zero, 


*  =  l*  u'  ~  V  -  ey  -  cos/f) 

~  ur  \  p0  -  e  j  ~  eA  smj3) 
r'*'  (?.’ey-e*  s/^] 

•  ur'^'  (i/c  ~eH  ~CA  ca3/f) 

e'A\(v'cosf-i-ur,sin,6). 


>v  taking  any  product  of  v  ,  *r  ,  <p  ,  ir' ,  cs' ,  <p' ,  v"  and  ur" equal  to 
it  follows  that 

i*'),-:ya  -<y]i 

=  f  ’*  -  (V ,x  *• 

~  Z  lr'  "e^  cos/;  -  2«r"(fre  -  e}  -eAsnyS) 


~  2  ir'{eA  cos, 5  -cA  0*  smti+e') 


(9) 


+  2  (e A  sm  S *  s.  d  cos  3  +  e'  ^ 


+  2  <p 


•y*“ey'ieV  *,n/s) 


-#o  ~^](Gy  r  e?  C0SA) 
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•  Let  dsa  and  ds  be  the  lengths  of  an  elementary  fiber  before  and 
after  deformation.*  Then, 


fe-J- 


and  the  longitudinal  (tensile)  strain  of  the  fiber  is 


ds  ~  ds 


'ds  y/cc  Sg 

dxj/  W  ■ 


T - - 


~ 2  <+<#;)'+&  f 


where  the  last  step  follows  because  ds/ds0  is  approximately  equal  to  unity. 
When  Equation  (9)  is  substituted  into  this  expression  and  it  is  assumed 
that (yj)* ■I’fya)*  <K  1  >  the  following  equation  for  e  is  obtained: 


The  end  points  of  an  elementary  fiber  in  the  two  adjacent 
sections  are  taken  to  have  the  same  if  and  the  same  C  • 


1 1 


"  -ey-*A  cos/3)  -  ~  sA  sm/3) 

+  ir‘  (e^  y-  cos /3  -  <cA  /g'  s/nfi ) 

+  ur  1  (e'j  t  e'A  s m  /3  +  eA  /S  'cos/3) 

+  ~  cy)(ef  +  e’t  *’”/*) 

-(fa-e)  *  e/  c*5/^ 

*{&“  eyf  +  ti*'erf}  /*']  • 

The  expression  given  below  for  6  is  obtained  by  letting 

q**t  +  n 

and  by  using  Equations  (1)  and  (2)  to  eliminate  y0  and  y0  . 
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6  =  u.'  -  (rj  -  eA  )  (ir"  COS +  ur"  smfi) 

-  £  ( ur ''  cos /S  -  v"  s/rt /. 3 )  +■  e y  ir'  +  e'^  ur  ' 


eA  (ir'COS/S  +  ur‘  5 m 

*  f'  iq  fe  j  cos/-  e'y 


/$)■*-;&  'cA  (ur  'cos /$-v  s/nfi)  ^ 1 
sm/3)- £  (e'y  cos/3  -he'^  s/n^-Ae^) 


The  following  elastic  properties  of  the  section  will  be  used  in 
the  subsequent  development: 


Figure  2.  CROSS-SECTION  SYMMETRY. 

K*Ta 

A 

Elt  =  f  C,X  E  d A 

•'A 

EI2  z  fjt]-eA)3  EctA  =  jf  q  (q  -  eA  )  E  dA 


EB^f  q(q‘+t‘-**)E<eA 

JA 

-  jT  ft**  Z'U*!  ~eA) £eCA 
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ES,  *  f  -  £x'  t*  )  E  ctA 


G  J  =  conventional  torsional  rigidity 

The  total  tension  load  applied  to  the  blade  section  in  the  X- 
direetion  is  denoted  by  7~  and  is  given  by 

T  =  f  £  6  cLA  . 

The  substitution  of  e  from  Equation  (10)  into  this  expression  yields 
u  =  -  c'v  ir  -  ur  '  -  e'A  (ir' cos/?  +  ur’  s/n j$) 

-  f'eA  {u.-  loss'-  r  t>m^ 

-  4>  [e'y  cos ,e-  ey  sm/S )  +  #A  ft' 

Where  the  average  tensile  strain  is  defined  by  e~~  £A  •  This  expression 
can  be  used  to  eliminate  a/  from  Equation  (10),  giving 


e  -  t  y  COS/3  r 


-  C  l.i’-  ccs  3  -  ir  '  s  m 

(ID 

r* 

*  °  |  (’I  ~  eA smfi) 

-t,  (ey  COS/3  *  e  ^  ±,n  f  *■  e  t)  +  £  z-  £*)/*'}  . 

Let  M,  be  the  flapwise  bending  moment,  M2  be  the  chordwise  bending 
moment,  and  <5  be  the  torque,  all  referred  to  axes  passing  through  the 
elastic  center  of  the  section.  M,  is  positive  when  it  tends  to  compress 
the  upper  surface,  is  positive  when  it  tends  to  compress  the  leading 
edge,  and  Q  is  positive  when  the  torque  applied  to  the  inboard  side  of 
the  section  tends  to  make  the  leading  edge  move  up.  The  bending 
moments  are  found  by  multiplying  the  stress  by  appropriate  moment 
arms  and  integrating  over  the  section,  giving 

=  -  f  C  E  e  d  A  < 1  2) 
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=  -/  if  £  e  dA 


It  may  be  noted  that  the  tensile  stress  £e  of  a  fiber  has  components  equal 
to 

-  EeZ(/B'+f) 


ana 

£  e  ‘  (*?  /V  f 

acting  in  the  rj-  and£-  directions,  respectively,  on  the  inboard  portion 
of  the  blade  at  a  given  section.  The  torsional  moment  due  to  these 
components  must  be  added  to  the  St.  Venant  type  torsional  term  &  J  <p', 
giving 

Q  =  GJ<f>'  +  f  ZE^tft'+t')  oCA 

(14) 

t ^  rfE  e  \(t\-zt)/S' +  rf<f>' ■  dA 

By  substituting  Equation  (11)  into  Equations  (12),  (13),  and  (14)  and  using 
the  sectional  elastic  properties  defined  in  the  last  paragraph,  the 
following  results  are  obtained: 

=  EIf  |  ur  "  cos  ft  -  ir"  sm ft 

(15) 

■t-  4>' (s'y  cos  ft  +  e  '  s/n  ft  +e't)  • 

Mz=  EIZ  iyr"  cos  ft  -  «r"s>n  /5  -  cosjS-e'yStnft)  -  (l6) 


-  TeA-E62ftft' 
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(17) 


<?  -  |G  J  +  *  (e'  cos/3  -  e'y  s/nft)(E  B2  -  e{  E I3) /S 

'(E6,-et  EBt)(/'!‘ ]•*'-  eA)f' 

-(EBj-  et  Elftft'  (ts"  cos  ft*- «r  "  stn ft)  . 

Let  dy  and  dy  be  the  y-  and  y-  coordinates  of  the  local  axis 
of  built-in  thrust  of  the  undeformed  blade.  Then, 

^ y  —  d  y  —  e  ^  cos  ft  and  6  y  -  ^ y  —  e  g  s  t  n ft  . 


Hence, 

e'y  cos ft  *■  e' y  smft*e'^-  dy  cos  ft  +  d'y  s/nft 

e'y  cos  ft- e'y  sm  ft  -  d'y  cos/3-  dy  stn ft -  et  ft' 
and  Equations  (15),  (16),  and  (17)  may  be  written  as 

=  E  If  cos  ft  -  zr"  sm  ft  +  1>'(dy  cosft -td 

Mz  =  £ I2  [r"  cos  ft  +  us"s/nft  -  <p  '(d'y  cosft-d'y 


(18) 


(19) 


T&a  (E  5 2  ~et  Eft  )  ft  ft' 


and 
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q  =  |  GJ  +  T-&*  +  (dj  cos/3-d'y  stn^)(E  Bz  -  £  1 2)  /' 

+  (EBr2et  £B2  -  e/  £IZ)(/?')Z  j  <P  ' 

(20) 

+  T(£\-et  eA)/}' 


-(EBz-et  £  Iz  )  jS '( v"  cos/2  +  ur  “  s/n/2 )  . 


If  the  axis  of  built-in  twist  is  a  straight  line  parallel  to  the  blade  refer¬ 
ence  axis,  d'y-d^O  ,  and  several  terms  are  zero  in  the  above 
equations.  When  the  two  axes  are  also  coincident  (  et  *  0  ),  the  above 
expressions  reduce  to  those  given  by  Houbolt  and  Brooks  (Reference  3  ). 

MOTION  OF  COORDINATE  SYSTEMS 

In  the  present  analysis,  equations  of  motion  are  developed  for  a 
rotating  blade  by  finding  the  inertial  forces  acting  on  mass  elements  of 
the  blade  in  terms  of  the  motions  about  the  flapping,  lead-lag,  and 
pitching  axes  as  well  as  the  blade  elastic  deformations.  The  expressions 
for  blade  motions  caused  by  rotations  about  these  axes  are  very  complex, 
particularly  when  possible  offsets  of  the  various  axes  are  taken  into 
account.  The  required  expressions  for  these  motions  are  derived  in  this 
section  by  a  series  of  transformations,  each  showing  the  effect  of  a  parti¬ 
cular  rotation  and/or  a  translation  of  the  origin  of  the  axis  system. 

The  first  six  coordinate  systems  which  are  used  in  finding  all 
the  blade  motions  exclusive  of  elastic  deformations  are  shown  on  Figure 
3  .  Two  additional  transformations  will  be  described  later  which  give 
the  orientation  of  the  blade  reference  axes  relative  to  those  shown.  The 
°r  *x  2x  J'x  system  shown  on  the  figure  is  an  inertial  system,  fl  is  the 
constant  rotational  speed  of  the  rotor  about  the  j-r-or  ^f-axis,  in  radians 
per  second.  The  y1  system  rotates  with  the  rotor.  Let  U;  ,  Vi  ,  We 
be  the  linear  velocity  components  of  the  origin  0 ■  in  the  directions 

and p-,  <}~. ,  r{-  be  the  angular  velocity  components  of  the  0t  ici y£  j,.  system  in 
the  z-t  g.  , y.  directions  {i  -  1,  2,  3,  4,  or  5).  Then, 

Ui  •  vt  m  W,*  f,  *  ff*  °  , 
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The  02  £zyi}2  system  is  obtained  by  rotating  the  Ofz1y^y1  system 
about  the  negative  yz- axis  (which  is  parallel  to  the  yf-axis  and  located  at 
x\~  rF  )■  This  rotation  angle,  p  F  ,  is  the  flapping  angle  of  the  blade,  and 
it  can  be  seen  that 

u2  =  0,  V2  =  flrF  , 

Pz  =  A  =  ,  r2  =  II  cos  ^  . 

The  OjZjXjfj  system  is  obtained  by  rotating  the  Oz%xyzyJ system 
about  the  ft  axis  through  angie^,,  the  blade  lead  angle,  and  then  trans¬ 
lating  the  origin  in  the  xzy2~  plane  to  zz - rpD  and  y  --  ep  ,  so  that 

V3  =  Pz  casA  i‘fz  s">  A 

=  il  S/n/SF  cos  ^  -fa  stn  /Sp 

h  =  C05A  ~Pz  SinAx> 

=  -n  s//7/a  eos/*i> 

r3  r  '*  *  11  COSJF 

U3  =  uz  cos  +  V3  sm/S^  -  r2  (eM  cos  ^  -  rp2>  ) 

=  n^rF  s/n  /$D  -  (es  cos ^  - rFM  sm  ^  )  cos^  j 

V3  =■  V2  cos/j,  -  uz  smjSD  r  rz  (es  *  rfJ>  cos/5^ ) 

~  ^  {  tf  cosA  r  (es  s,n/^s  ^  cos.#j,)  cos>AP  ' 

. 

W3  ■'  Wt  ~Pzev~  rr-S  =  n  e2>  "Ar  ^  . 

The  0V  system  is  obtained  by  translating  the  origin  of 

the  03  x'}y}h  system  to  *,  =  rsp  ,  y3  --e^  and  fi:Cp  ,  so  that 
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pf  -  p3  -  a  5/n  0F  cos  SD  -  8p  Stn 

9-y  =  ^  5  S,n^D  ‘A  COS  A. 

rf  -  r,  =  .ft  cosjp  t-  sD 

Uv  z  US  *  ^  Ca>  eSP 

-  &{(/>  -  Cp  sm  ft F  )  s,nj3D  -(e  DP  +  e0  c os/*  -  ^  S'"/A°a} 

”  A  cos  A  ~/$d  eDP 
^  =  V3  ^  rj  ri>/5  _  Pi  C'p 

=  ft{fr>  -  Cp  s/«/p  cos/,  v  ^  *  e,  s,nfa  +  rrj)  cosjSD  )  cosfc  | 

+  p  CP  S,n/z>  */4  rDP 
Wt  r  *  f3eOp-?3rPP 

=  n.(e9  +  eop  cosfa  *  rsp  sin/$s)s,n /?p  *•/*  ^  *  V  CosSiTej>P  S//7A)* 


The  y5  system  is  obtained  by  rotating  the  0+  x+  yy  system 
about  the  £¥~axis  through  an  angle ftp  ,  the  blade  root  pitch  angle,  so 
that 
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U c  =  U9 


=  a  \(rF*rrD  cospF  -  Cp  sm fip)  sm/SD 
'(evP  +  e*  cos/Jcosa} 


*F  ^ P  COS/v~  A>  e3P 


V5  =  v;  cos/3p  +  w¥  sm  j3p 


=  n 


{fa  +  rFD  C05A“  °p  s/nA)  C0SA 

+  (V  +  ei) 5/ri A)  cosa}  cosA 

*  A  "■  ej>P  cosA  *  ^  Sin^x>)sm/SF  S,n/sp 

+ fa  •  ^  s//»A  cos/5p  +  (rpj>  +  rsp  cos/3 s  -  eDP  s,n/SD)  sm/3p 

+  fa  r*P  Cos/P 

n/5  =  Wy.  cos/p  -  V \  S  m  /Sp 


=  -O. 


(ej>  +  evp  cosA  +  rDP  s,n  A  )  S!n  fa  cos  fa 


fa  +  rFP  COS/f  '  CP  Sm/p)  cos  A 

-fa  r~~  s,n 


Z>  'DP  /A 


*  fap  +  eD  s/nA)c°5Aj  $,nfa 

_j 

*  fa  \fap  +  rDP  cosA -  s,nA)cosA  ~  CP  s,”/v * "»  a{ 
Vs  =  ^  +/p 


=  -°-s//7A  cosA  -  a 
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'h  '  i  \ 


os  o 


■Cl  (sm  3p  bin  Jv  <-cs  3p  -  cos  SF  S'n/3p  ) 


3p  CCS  3p  CCS  Op  S/n  3 P 


y\  ~  n  i-os  .i  -  a  6' n  O 

>  ^  p  r  u. 


,p 


"  n  (cos/A  CCS  3 p  -  Sin  y.  Sin  Sp  b'n/3p  ) 


'/f  cos,  ^  s  ">.  JP  ~A  C**,3P 


These  six  relations  mav  be  written  in  matrix  form  as 


where 


1 

'5. 

0 

4 

1 

/5 

5l 

r> 

^4- 

4 

0 

n 

c 
‘  s 

rx 

C 
‘  r 

0 

A 

= 

5 

~v 

u 

0 

< 

,  > 

4> 

s„ 

* 

r- 

‘  V 

0 

\v  s 

—  3  — 

IV 

3, 

0 

(21) 


*cs 


-  cr  -  *  P  ^  n  P  ^  ~ rp 


c-i  y.  yp  -  S/7  3p  c,m  3 D  S'”/3p 


>  ? 


2u  =  (rF  +  rFO  COSA  "  Cp  S/n/3F)S,nA  -  (efip  r  eD  COS/j>  )  CC,SA 

SV  =  jo>  -  rFD  cosA '  cp  s,nA  '  cos4  *  fo*  "•  **  y°sA-J  cos,A 

*  A  *  evP  cosA T  r»P  s,nA ) s,n A  s,nA 

V-j (tf  r  rFv  cosA  ~  cp  s'nA)  cosA  +  v "  eD  s,nA )  C0SA  } 

*  A  "  ei>p  C0*A T rDP  5,nA)  s,nA  cos/*p 

F?  -  -  S/nA 

Ff  a  -  «sA  cos^P 
Fr  =  CosA  S,n/3P 

Fa  =  ~Cp  cos 

FV  '  (rFD  +  rvp  COS/p  '  ePP  iinf» )  s'”/p  *  CP  S'” A  cosA 
Fw  =  (rFD  +  r2>P  C05A  '  ePP  5/0 A  )  C05  ,SP  "  CP  SmA  5,n/SP 


-  btn  Ap 

Dr  r-°’-Sp 

II 

rf 

A  "  e2>P 

Sr  V  C0SA 

~  ~  rPp 
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It  follows  from  the  preceding  matrix,  e.  g.  ,  noting  D  -- d ,  that 


As  mentioned  previously,  the  1'inal  equations  of  motion  which 
are  obtained  tor  tinding  natural  vibration  frequencies  and  mode  shapes 
ot  rotating  helicopter  blades  are  simplified  by  the  omission  of  terms  of 
negligible  importance.  These  equations  could  be  obtained  by  first 
working  out  exact  equations  and  then  linearizing  them.  However,  a  con¬ 
siderable  simplification  can  be  obtained  in  the  analysis  if  terms  are 
dropped  in  the  course  of  the  derivation  when  it  beco.mes  clear  that  they 
will  lead  to  negligible  terms  in  the  final  equations. 

The  first  time  derivatives  of  the  blade  flapping,  lead-lag,  and 
pitching  angles  (1.  e.  ,  ,  /SD  ,  and )  appear  explicitly  on  the  right-hand 

sides  of  Equations  (21)  and  (22)  and  also  enter  implicitly  through  the  time 
derivatives,  S?  .  ■“*  ,  ,  Zy  ,  etc.  The  ratios  of  these  time  derivatives 

to  the  rotor  angular  velocity  are  assumed  to  be  quantities  of  first  order 
which  are  small  compared  to  one.  This  assumption  might  be  justified, 
tor  example,  by  assuming  sinusoidal  oscillation  of  the  blade  at  a 
trequency  of  lOH  which  involved  a  blade  root  flapping  motion  of 

-  0.  5  deg.  =  0.  00872  rad.  In  this  case,  Sp/a  -  (10A  )(0.  0087)//!  =  .  087. 
Lower  values  of  this  ratio  would  be  obtained  at  lower  oscillation  fre¬ 
quencies  and  amplitudes  oi  vibration. 

It  Equations  (21)  and  (22)  were  used  in  deriving  expressions 
tor  the  absolute  acco’ orations  ot  a  n  ass  point  m  a  blade,  a  number  of 
terms  would  be  obtained  involving  products  of  <5r  ,  ,  and  .  These 

products  could  be  neglected  because  they  would  give  higher-order  terms 
in  flf/n  ,  fiD  A  .and  6p,i \  it  the  expressions  were  written  in  a  suitable 
nondmiensional  form.  Some  additional  terms  would  be  tound  in  the 
expressions  for  the  absolute  accelerations  in  wfhich^  ,  ,  and/^,  would 

enter  linearly.  These  terms  could  be  identified  as  Coriolis  accelerations. 

The  Coriolis  torces  acting  on  a  rotating  blade  produce 
couplings  between  vibration  modes  which  are  uncoupled  when  the  blade 
is  not  rotating.  These  intermodal  couplings  are  unusual  in  that  they  are 
oroportional  to  the  time  rate  of  change  of  the  variables  ( ,  and  ) 
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and  they  are  unsymmetrical.  An  interesting  consequence  of  couplings  of 
this  type  is  that  the  exact  solutions  for  the  "natural  vibration  modes"  of 
the  rotating  system  would  give  motion  in  which  all  parts  of  the  elastic 
body  would  not  move  exactly  in  or  out  of  phase.  However,  these  effects 
are  small  in  the  case  of  rotor  blades  and  can  be  safely  neglected  when 
computing  their  vibration  frequencies  and  mode  shapes. 


It  follows  that  the  terms  involving  fa  ,  ,  and  f>p  in  Equations 

(21)  and  (22)  will  have  a  negligible  effect  in  the  final  equations  and  can  be 
omitted  in  the  following  derivation.  Thus,  Equations  (21)  and  (22)  are 
simplified  as  shown  below: 


■Ps 

r 

5? 

h 

54 

r5 

< 

>  II  < 

> 

US 

Vs 

5, 

w5 

S* 

*  - 

V 

'fir 

fs 

Ff  ftp  *  D 

+5 

+  ^rPn 

<  . 

>  ~  < 

..  f 

U5 

F V  Pf  +  Pv 

i 

Fv/>F  +  ^vPd 

to 

FW  Pf  +  Pv 

to.  -* 
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(24) 


The  (  %, ,  jf  1 1  yi) >  » Ux  >  3  > !/j  >  + » y* >  y? ) »  and 

/  v  >  i  t.  \  p  vaci  cvofomo  U  a  v  n  on  inf  nrl  ir»  f-1*  « 

\  i  Jg  t  0$t  ^  ***•  *  «  ww**  wa  uvi  vicxxaxc  lac 

linear  and  angular  motions  of  the  blade  due  to  shaft  rotation  and  motion 
about  the  flapping,  lead-lag,  and  pitching  axes.  The  zyf  axes  system 
which  is  used  as  a  reference  system  for  the  blade  has  a  nxed  orientation 
relative  to  the  Z9ysy5  axes  system,  as  indicated  on  Figure  4. 
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The  zyy  system  (with  the  2T-axis  being  the  blade  reference 
axis)  is  obtained  by  rotating  the  Zsys%s  system  first  about  the  ^5-axis 
through  an  azimuth  angle  dA  and  then  about  the  negative  y-  axis  through 
an  elevation  angle  Qg.  Thus,  the  rotational  velocity  components  of  the 
Zyf  system  in  the  z-t  y-  ,  and directions  are: 

P  *  (Ps  cos  9A  +  ys  Sm  &A )  cos  &£  *  sm  0£ 

9  ‘  fs  C0S  eA  ~  PS  S' »  9A 

t  -  r$  cos  9£  -(^>g  COS  Ba  f  fs  s/n  &A)  s/n  &£ 

They  may  be  written  as  follows  using  Equation  (23) : 


where [t]  is  a  3-x-3  matrix  with  the  following  elements: 

=  cos  &A  cos  e£  rxz  =  sm  eA  cos  e£  rXi=  s  m  e. 

r2t  =  -  sm  eA  cos  9 a  (rzs  =  0) 


fJ}  -  -  cos  &A  s/n  Q£  r32  -  -  Sin  &A  sm0£  T}3  =  cos  &£ 

Similarly,  the  translational  velocity  components  of  the  origin  of  thexyy 
system  in  the  X-  ,  </-  and  directions  are  given  by: 
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The  xj  y  components  of  the  angular  accelerations  of  the  x.<jy 
system  and  the  linear  acceleration  components  of  the  origin  of  the 
/jv  system  can  be  obtained  from  the  z,  u.  components  by  the  same 
transformation,  giving  * 


The  total  flapping  angle  and  flapping  acceleration  can  be 

written  as 

fip  ~  ^r^osal.  r  ftp 

P  f)oscil. 

where  /6F  is  the  steady  or  average  flapping  angle  and  (/^p)wlU  is  the 
oscillatory  component  of  the  flapping  angle.  A  similar  breakdown  into 
steady  and  oscillatory  components  can  be  made  of  jSB  and  p  and  the 
elastic  deformations.  It  is  required  that  final  vibration  equations  be 
linear  in  the  oscillatory  variables.  This  will  always  be  a  reasonable 
approximation  for  sufficiently  small  oscillations. 

In  addition,  the  mean  angles  j$F  ,  fiD  ,  and^ ,  the  azimuth 
angle  and  the  elevation  angle  (Ge)  are  considered  to  be  small  which 

leads  to  simplifications  in  the  equations  for  the  acceleration  components 
obtained  in  the  next  section.  Although  these  constant  quantities  are  not 
considered  to  be  as  small  as  the  oscillatory  quantities,  they  are  assumed 
to  be  small  enough  so  that  their  squares  and  products  can  be  neglected 
when  they  appear  in  the  coefficients  of  the  oscillatory  variables.  When 
the  expressions  for  -p  ,  ^  ,  r  ,  etc.,  and  simplified  accordingly,  for  use 
in  the  next  section,  they  take  the  following  form: 
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p  =  PL  (d£  fpF ) 


l  --  A  A 
r  --  a 

I 

t 

U  =  n.(r9^i-rp6^-ep) 

v  s  n  (ep  A  *  r/»  r  eo  A  ‘  Cp/3f) 

W=  Pi  p  A  *  ep  A"  "  r  p  ftp ) 

>  (25) 

f>  =  "A  A  A 

?  --  -A  +AA'  A  A 

r  =  AA^A'A  A  i 

M  ~pF  [Cp-(rp  -  rp)  A  '  +  A>  ("  **/>  *  DP  6 A  ) 

v  =  a  {a A  A, My  a ) A  ■ *  A,  (y  A  *  y) 
w  =  A  (A  A  *  rp~rp  -  A*  4) +  A  (e»  A  “  A  J 


The  blade  is  a  beam  in  the  Zyf  system  cantilevered  at  K  =  0. 
Consequently,  a  mass  point  in  it  can  move  relative  to  the  zy  y  system  as 
described  in  the  section  on  the  elastic  behavior  of  a  twisted  beam.  This 
motion  will  be  taken  into  account  in  addition  to  the  velocities  and  accel¬ 
erations  given  by  Equation  (25)  in  finding  the  inertia  forces  acting  on  a 
mass  point  in  the  blade. 
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In  the  present  case  the^y^-  axes  are  the  reference  system  for 
the  blade.  In  the  section  on  the  motion  of  coordinate  systems,  expres¬ 
sions  were  derived  lor  the  aryy  components  of  the  linear  \olocity  of  the 
origin  of  the  xyy  svstem  and  lor  the  components  ol  the  angular  velocity  oi 
the  x.yy  system.  Expressions  were  also  obtained  lor  the  time  derivatives 
of  these  components.  The  approximate  lorms  ot  these  expressions  for 

p . W  .  given  in  Equation  (25)  are  to  be  substituted  into  Equation  (27i. 

Equations  3).  (4).  and  (5)  derived  in  the  section  on  the  elastu  behavior 
of  a  twisted  blade  express  the  coordinates  ol  a  mass  point  oi  the  blade  m 
terms  ol  the  elastic  deformations.  These  expression-  :or  i  .  j  .  and  } 
are  also  to  be  substituted  into  Equation  (27). 

The  underlined  terms  in  Equation  (27)  represent  Coriolis 
accelerations  associated  with  the  velocities  oi  a  mass  point  relative  to 
the  zyy  system  arising  tram  elastic  deformations  ol  the  blade.  As 
discussed  earlier,  the  inertia  couplings  produced  by  Corioli-  lore  as  hav«- 
only  a  small  influence  on  the  vibration  trenuencies  and  mode  shines  ol 
rotating  helicopter  blades  and  can  be  neglected. 

Some  new  notation  and  some  changes  in  r.otatioi  v  ill  ).«  made 
in  carrying  out  the  above-mentioned  substitutions.  The  s\  mboi  r  «  jii  no 
longer  be  used  to  represent  the  omponent  ni  a*  gular  vein*.  tt%  but  will 
be  denned  b\ 


»"  =  "  *  xe 

The  physical  sigmticance  oi  the  redetiried  ouantitv ,  *  an  b*-  seen  most 

readilv  lor  the  case  where  40  <?4  o  where  1 1  is  the  rad»ai 

distance  ironi  the  rotor  shait  to  a  blade  »e.  t . > >n  (who  t.  s  »  her,  r  •.  rn  al 
to  the  blade  reierence  axis). 


The  deiiectionx  o’  t he  eia x 1 1 <  «  • ,  ■  ,n  Mif  ,  u  11  }  '.ret’  i ,  > , 

which  ar»  produced  b\  the  total  ^umoireo  i,  ..u  .  er-.ij;  o  "  ’  r-  m 

are  denoted  resprv  t,v  t  1\  In 


It  t hen  tedb ,i»  .  :  •  f. i • ,  E •  (  i-' T ■  ■  -  1  4 1  <r.d  oi  n.  ii 


j  t 


<2h! 
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*  t  . 
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Only  very  small  oscillatory  motions  ot  the  blade  about  the 
blade  steady  contiguration  are  considered.  Let  irr  ,  $  ,  sf  ,  aD  and/<<, 

in-  the  steads  portions  ot  ly.-ey,  <t>  ,  "<.,^.11111^.  respectivels .  while  the 
latter  set  ot  s\  minds  w  ill  represent  trnm  now  on  anl^  their  oscillatory 
portions.  (  The  subsi  ript  is  omitted  troni  the  symbol  Asse//  used  in 
the  previous  ms  turn  leaving  the  symbol  tor  the  oscillatory  tlapping 
allele  ) 


Bs  malms;  tn**  substitutions  svbi»  h  have  iieen  disc  u  *<  set!  on  the 
n^ht-hand  sid«  *1  Equation  t  J!  7 1  observing  that  *  .  .  •/,  <>,  and 

dropping  any  protju*  t  ot  the  <ui  illatorv  displa.  or  their  time 

der  iv  itivi'..  ftn*  !oilo\s  results  are  obtained 
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Here  the  total  components  at  the  absolute  acceleration  are 
denoted  by  :it  ,  <ty  .and  «j,  ;  the  oscillatory  components  by  **  .  <l  ?  ,  and 
4  y  .  and  the  steady  components  due  to  A  by  .  and  ay. 

As  discussed  earlier,  <?4  and  and  the  steady  portions  at  the 
nandimonsional  variables  are  assumed  to  be  small  quantities  although 
not  necessarily  as  small  as  the  oscillatory  variables.  The  squares  or 
products  of  these  quantities  are  neglected  in  the  coofiicienta  oi  the 
oscillatory  quantities.  However,  products  oi  these  quantities  have  been 
retained  in  the  nonoscillatory  parts  of  the  expressions  which  enter  into 
the  determination  of  the  average  deformations. 
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It  can  be  readily  seen  that  the  inertia  forces  produced  by  the 
<  ,  component  ot  acceleration  are  essentially  in  the  axial  direction  and 
are,  m  general,  multiplied  by  oscillatory  angles  or  moment  arms  in 
the  \ibration  equations.  Hence,  after  inspecting  Equation  (30),  it 
appears  reasonable  to  approximate  by  the  familiar  expression 

^  (33) 

dropping  terms  alroad\  multiplied  by  oscillatory  quantities  plus  a  few 
additional  small  terms.  These  additional  terms  which  are  dropped 
would  be  much  smaller  than  the  one  given  above  except  possibly  those 
very  close  to  t»(e  blade  root. 

A  similar  simplification  of  the  expressions  for  and  E is 
not  possible  since  no  single  term  predominates  for  these  components 
and  the  corresponding  inertia  forces  act  normal  to  the  blade  axis. 
However,  when  the  above  formulas  are  used  in  obtaining  the  blade 
inertia  loading  later  m  the  analysis,  all  the  underlined  terms  in 
Equations  ( <  1 1  and  (id)  will  be  dropped  by  considering  that  the  geometric 
eitects  ot  the  steady  delormations  ib  ,  <e-  ,  <P  and  ft  on  the  oscillatory 
motions  are  either  negligible  or  would  be  included  in  ey,  e y  ,  ji  ,  9A  and 
9e  ,  (That  i>.,  the  initial  blade  coniiguration  would  be  defined  including 
the  steady  deformations  produced  by  centrifugal  forces.) 


BLADE  INERTIA  LOADING 


The  inertia  torce  on  an  lntimtesimal  element  of  the  blade  is 
obtained  t>\  simple  multiplying  its  mass  times  the  negative  of  its  abso¬ 
lute  ,u  t  deration  as  given  In  the  expressions  in  the  preceding  section, 
fur  omv  I'liu-m  r  in  the  tollowmg  analysis,  these  infinitesimal  forces  are 
ml*  grated  o\er  a  section  to  obtain  expressions  tor  the  inertia  forces  and 
moments  ,s>  ting  on  the  blade  per  unit  span.  A  blade  section  is  con¬ 
sidered  whuh  perpendicular  to  the  t-axis  of  the  blade  reference 
y  stem  and  located  at  »  x^  in  the  unloaded  condition.  The  components 
<>!  in*  rtia  tori  es  and  moments  per  unit  span  acting  at  this  blade  section 
are  u r si  determined  yyith  respect  to  axes  yvhich  are  parallel  to  the 
'.lade  reteretic  e  system  but  yyhose  origin  has  been  displaced  a  distance 
t  1 1  ilong  the  r-u\is.  The  components  of  the  inertia  lorces  are  then 
in  axl  relatn  t«>  axes  which  are  parallel  and  perpendicular  to  the  blade 
c  hot .  i  tor  the  giyeii  section. 

I  bus,  with  c  being  the  mass  density  of  the  blade  material  at 
an  .  mmi  ni.iri  .  rus#i-st*i  tmnal  area  dA .  and  ait  Hu  .  a.,  being  the  absolute 
at  celt  rations  at  the  six  torce  and  moment  components  per  unit 

span  .i  rt  ; 


-  [«■*  p 

-l s* ' 


cCA 


d  A 


Py  -•  ~  f  Zt  *  *A 

f/zy}  '*> 
f  sxy 


dA 


(This  moment  compo¬ 
nent  is  in  the  negative 
y  -  direction.  ) 


The  moment  components  about  the  elastic  axes  o£  any  blade 
section  are  readily  computed  from  these  inertia  force  and  moment 
loadings. 


The  expressions  lor  <£x,  Sy  ,  and  a j, given  by  Equations  (31), 

(32),  and  (33)  and  the  expressions  for  <j  and  y  given  by  Equations  (28) 
and  (29)  are  to  be  substituted  into  the  above  equations,  and  higher-order 
terms  in  th<  oscillatory  variables  are  to  be  dropped  as  in  the  preceding 
section.  The  initial  position  of  a  blade  element  which  appears  in  these 
expressions  is  then  written  in  terms  of  the  section  coordinates  indicated 
on  Figure  1  as 


ifo  ~  e y  *  *1  cos /  '  S’ 

y0  =  ey  r  1  b,n/5  r  S 

It  is  convenient  to  define  the  following  sectional  inertia  pro¬ 
perties  for  use  in  the  further  development: 
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^  5  4r  f  ft  p  d  A 

r 

j  /?  dA  -  w  (e^  +  rjc  cos  jd) 

f  }o  Pd A  ;  **  (ey  r  % 


A'zik3  -  4* 


tvKx  =  L  j* ^0  dA  =  *n  (e^  -r  Ze y  rj c  smfi *•  ^  s*/')  /  *■  cos  /$) 

4 

tnKz  =.  J  <f*  /><*A  =  r?  (e^  -  -Ze^  /^c  cos/  ■>■  &]  cos'  ft  *■  smft) 
>r  H  ■U»  0  dA  =  »?  |e^  e^,  *•  (<2j,  s/^/’  r  e  cot, ft)  rjc 

+  (k2  -  A*)  Smft  COSft^ 

K  —  K {  *  K2  -  ^  4  *  ft  y-  e  *  *  ^  r]c  (ey  cos ft  ey  5,nft ) 


K,  '  K,  a  (**  -  */  )  cos2ft+  e’  -  e*  -2 rjc  (ey  cos/  -  smft)  . 

Also,  since  the  blade  is  assumed  to  be  symmetric  about  the  chord  line,  it 
follows  that 

j'  C p  dA  -  0 ,  J~  f]  C  p  d  A  =  0,  etc . 

Using  these  properties,  the  sectional  properties  defined  above  and  the 
expressions  for  ,  Zy  ,  and  5/  (excluding  the  underlined  terms  in 
Equations  (31)  and  (32)),  the  inertia  loading  components  become 


p„  -  fl  y*i r 


-  **  fa  +/D)(CP  +  ^  7c  s  /n/ )  Ar 

*-  0,  ^  fe;  *  £ 

-  *  js£  (r-rp)+  ey  *  yc  s/n/f  j/^ 

*  ^  [/•  -(ic  i»/)*-i4-{e?ts"/'(4'/>-''i 

+fl’w  (■=/■  -  /,  r/=  -/£  'i  *•  «y  *  7C  c»s/;  (35> 

-  7>t  -7*1  (ffc  COS/S  )<f> 

-  -m  •  r  -  rF  +  d£  Cp  ~/D  eMp  -  (0A  +/D)  fe  ^  *■  qc  cos/;J  fi '  p 

'lr'eE  (fiBp+  +  yc  cos/?; 

-  «N 

->f  s  ^  (r-rp)  +*;,  +qc  cos/  /p 
^  J 

"  A*  *  j  /y£  *  1-ay  +r]c  cos/)  /\ 


.?  -  »•  t>  -r.  i,  nS'ii  ->’t{ay+r[cccs/S)*'c-»!^’"lc(e',tos/a~<:rsin/?)^? 


-  »■  \T  i  j  -.  3"r <^p)  ie^  t  rjc  cos/3  )-(0A  +/SD  ) jtf  rCP(ey  r1c  s/,1Aj  ftp 


1  (r  ^  'i°ppXeyr1c  *  epp(e3^c  C0S/')J  A 


-  /y  ^ncL0S/3)*  0£  {Cy  +  qc  s/n/'j*-  K  /6p 

(9£  -  Sp)  r  y  e y  -r  rjc  s  ,n/3 j 
-nV  -  0,  rp  -  3}  rD  r-  *■  cos  ,£)  <^r 

-  a„,  |(ep  -  0,  rp  *,)  ?tcos£  //N  t')cos2/S+  ?c(ey  cos/S-e^sm, 

-  Cl  hr  rKey  -  /;c  cos/  ’  /.  -  ^  A  *  qc  smd)/3D 

"V 

+  j<2  -  *,  *■  *  *  A  '  ft  c<?s  S '  ’  ■ A  +$F  )(r‘rP  )(A  *  1'  s,n A  •  A 

-  Pl»\ .  -ft  r,-4  r^)  ve^  *  ?c  »  ^>(^^)r(ey^c«i/J,)| 


-  -  O.  rrj  t~  -ft  t,£,s/'  ■>}  -  Pirn  r  (ey  *■  rjc  sinfi) 


■j  -  -  Pi  »>r |  rt  -  xnc  s  ’/  ,'  S*j  "  AVr  {fty  -  ft  cos/  ' 
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In  the  expressions  for  each  of  the  components,  the  last  term  does  not 
depend  on  any  oscillatory  displacement  or  its  time  derivative.  There¬ 
fore,  these  are  steady  terms  and  will  be  designated  px  ,  ,  p,  ,  f x  , 

,  and  if.  for  easy  reference.  {Note  that^,=  px  =A’»?r.) 

In  the  subsequent  application  of  the  associated  matrix  method 
for  finding  vibration  characteristics,  the  blade  is  divided  into  spanwise 
segments.  The  treatment  of  individual  segments  employs  axes  which 
are  parallel  and  perpendicular  to  the  blade  chord  of  the  segment  (in  the 
unloaded  condition).  The  figure  below  indicates  the  resolution  of  the 
loadings  and  deformations  into  components  which  are  parallel  and  per¬ 
pendicular  to  the  chord  at  a  particular  segment. 


h>s> 

n 


Figure  5.  RESOLUTIONS  OF  LOADS  AND  DISPLACEMENTS  PARALLEL 
AND  PERPENDICULAR  TO  LOCAL  CHORD. 


The  deformations  of  the  elastic  axis  due  to  bending  are  denoted  by  <5, 
and  Sz  and  the  components  of  the  total  deformation  due  to  the  combina¬ 
tion  of  bending  and  torsion  by  8lt  and  82t  .  The  initial  coordinates  of 
the  elastic  center  and  the  axis  of  built-in  twist  are  denoted  by  (e,  ,  e2  ) 
and  (  di  ,  )  respectively.  These  components  and  their  derivatives 

with  respect  to  x  can  be  found  from  the  components  relative  to  the 
X'  ,  y-  ,  t-axes  by  the  following  transformations; 
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65,t  ’  5,  *  S',  >  or  O  ~  fat  »  M/*>  ur>  or  cos/ 

-  fat  ,  v,  v'  or  y")sm/ 
(sae  ,  S7  ,  s'  or  s")  -  fat  ,  ur,  ur'  or  ur" )  sm/ 

+■  (irt  ,  ir',  or  v")cosjS 
(pf  or  df)~  (e-y  or  dy  )  cos/  -  (e^  or  dy  )  s/n  / 

(ej  or  d*  )-  ( e y  or  d'y  )sn?/  -h  (e y  or  dy)  cos/ 


>  (40) 


J 


Expressions  given  below  for  the  components  of  the  force  and 
moment  'oadings  with  respect  to  axes  which  are  parallel  and  perpendicu¬ 
lar  to  the  blade  chord  are  obtained  by  similar  rotation  transformations. 

It  should  be  noted  on  the  diagram  that  jEj,  and  :ire  the  inertia  force 
loading  components  perpendicular  and  parallel  uj  ihe  chord  respectively 
but  that  ^  and  fh  are  the  moment  loading  components  taken  respec¬ 
tively  about  axes  which  are  parallel  and  perpendicular  to  the  chord. 


s  Sy  cos/  -  sin / 

-  -  ”1  S  -  nj  n  <A 

1C  I  c  ' 

-  *  \(r -  %  *  &£  V/4  cos/~  [dA  tA  ^  5//7/  *•*  +  %)}*? 

*  ^  '  5/V“  e9P  (e£  cos/  -  eA  sin/)  -  0£  (et  e  rjc  )  j  /£ 

"  *  jlf'* -/ )(*, r  *>»/  *  ^  COS/  )rc2tljc  J/p 

-fi  n^5J^os/  •> -r  S  -  S1t  s//7 1  £  -  (r]c  s/s?  */ / J* 

-  A*  r  vcos  J  )  -  fL*  Sm/ '  /5„ 

-A  nr  cpcos /  qc)  cos  2 /i  -  e ,  sm  2/ -{6£t/F){r- r/ s/n/ ~ jyT, 

r-toyCos/  '  p  o>n/  (4!) 


*10 


fi 


1 


-  -  rrj 


Py  r  Py  C°%/ 

*2t  '  *  *  $£  CP  'A  £*P]  *"»/  +  (94  V*  )(CP  C0£/S'h  *'  '}  ^ 

-*?\(r-rD)  cos /  +e-j,f>(Q£s m/  *  6A  cos/)  ~6£e1  /D 

+  *t^r-rp)(e£cos/-&A  sm/)-hei J /, 

+fl  *  7*  ^8}t  cos3/  -  Ste  sm/  cos/  -  (qc  sm/  cos /^  <)  ' 

-fl  T^f  (s/J)/.  r D  cos/)/D 

-R1  *1  | epsm/  +  +  )  sm  2/  +  ef  cos  2/  ^(^//(r-r/  cos/S^ 

+  fy  sin/  +  ?y  cos/  ( 41  ) 


?-y  COS/*  ~  fy  *"*/? 


-fL  m  r  ^  cos/  -  £  sm/ 


(43> 


h  * 


C„  5 


h  S>",^s  +  h  COS/S 

=  -fllntr  Szf  *■  fy  sm/  ¥■  f^  cos/. 

When  the  quantities  Cy  and  are  defined  as  below, 
ey  +  1c  cos/8  =  ( e2  r  1 1  )  c°s/  ~  ef  sm/ 


(44) 


Cy  -  ey  +  rjc  sm/  =  (?xrqc  )  sm/  +  cos /  t 
it  follows  that 

^  cos/  -  Cy  sm/  =  (e2  +qe)  cos  2/  -e/  sm  2/ 
\^y  itr>/  r  Cy  cos/  -  (ea  *•  qc )  sm 2/  *■  e{  cos  2/ 
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j  -  -  »’  ^  "  tje  '  &1t  r  "T  C  2  S2r  f).  ez  )  $ 

-»>{(r~r^+e.  Cp  -4  Cj>p)  Cy  ~  (9a  *  /D)(X  -  Cp  Cy  )}  /; 

*  ™{vr  -4  *  dA  e  dp' *  eBP  Cy  ]}/z 

-  '  (t)A  Cy  +  ee  Cy )  +  /cj-  ftp 

-  nl  rrj  -|(e;  trjc)  cos  2 -  ef  sm  2/S 

*  (eP  *  &A  rp  -/4  rD  )  cos/S -(6 £  rjSp)  r  stn/5 

-  \l  /»  j(e3  *■  qc  )  sm  2/S  v-  e1  cos  2/S 

*  tp’6A  'p-A  rj>)  S'"/3  *(6t  *fip)  r  C0S/^}  6Z 
~  ft  nr  |(c^  -  &A  rp  -£  rp  )  ^  COS /S  *(*]  -  */  )  cos  2/S 

r  rjc  (e ,  cos  2/S  -  c,  sm  2,4 )j  <P 

-  £i  mr  Cy/3F  ■*■  fi* m  rD  /SB 

-±l\  |<2  -K,  +ep£y-  fe£+/p)(r- rp  ’  Cy\  fp+h 


(45) 


42 


j  i ,  %  -  •  f  i  *'  i  ,i  *  l  •  * ' ,  *  ‘ '  r  ttit'iti-'  1 1  *•  i 


II 

-  fl1  W  •" 

/  X 

11 

(ep  -  dA  rp  - 

II 

-  11  77?  ($£  T-  y3p  )  i 

•ya> 

* 

it 

~fl 2  777  | ^  +  (ep~&A 

II 

2 

-fl  77?  r  Cy 

ii 

-fltnr 

(46) 


In  the  computational  setup,  the  blade  is  divided  into  a 
number  of  short  bays  and  the  mass  of  each  whole  bay  is  considered 
to  be  concentrated  at  the  mid-span  of  the  bay.  By  considering  th  to  be 
the  mass  of  one  bay  (equal  to  the  product  of  the  span  ot  the  bay  and  the 
average  mass  per  unit  span  of  the  bay),  all  the  above  expressions  for 
the  p's  and  p  s  express  correctly  the  inertia  loads  ot  one  bay. 


ACTING  BENDING  AND  TORSIONAL  MOMENTS 

In  this  section,  expressions  are  derived  for  the  bending  and 
torsional  moments  acting  at  any  blade  section  within  a  given  spanwise 
segment  or  bay  of  the  blade  as  a  lunction  oi  the  loads  applied  at  the 
outboard  end  of  the  bay  and  the  deformations  in  the  bay.  These 
expressions  tor  the  acting  moments  are  combined  with  expressions 
for  the  elastic  restraining  moments  in  a  later  section  in  order  to  obtain 
differential  equations  from  which  the  total  deformations  across  each 
bay  can  be  determined.  This  solution  tor  the  delormations  in  an  "elastic 
bay"  is  carried  out  considering  the  bay  weightless  because  the  effects 
of  inertia  forces  are  introduced  as  concentrated  loads  at  the  junctions 
between  bays. 

(  1 1  1  t 

Consider  a  short  bay  of  a  blade  between  »•  -  r'  and  «"  _  *•' »r 
with  constant  /S  ,  arid  ey.  These  constant  geometric  parameters  are 
designated  Ai'  ,  ey{t‘  ,  and  e'y  .  The  total  load _acting  on  the  outboard  end, 
where  r  =  r'>’  ,  has  the  force  components  Px  ,  Py  ,  and  Py  and  the  moment 
components  ,  and  .  These  six  components  can  be  found 
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respectively  b>y  integrating  the  inertia  loadings  pt,  p3  ,  p  ,  j^.and^ 
irom  r-r‘*!  to  the  tip.  The  pitching  moment  is  taken  about  the  intersec¬ 
tion  of  the*  -axis  and  the  plane  of  the  outboard  end  of  the  cross  section. 
(As  usual,  Qy  is  positive  in  the  negative  y  -direction  while  the  force 
components  and  other  moment  components  are  positive  in  the  positive 
directions  indicated  by  their  subscripts. )  When  all  terms  containing 
any  oscillatory  displacement  or  its  time  derivative  as  a  factor  are^ 
dropped  out,  the  load  components  are  designated  by  4  >  4  >  py  >  an^ 

and  may  be  called  the  steady  portions  of  ,  Py  ,  ,  Qx ,  ,  and  , 

respectively. 


A  local  coordinate  system  (  z(r),  y(r )  ,  y{r)  )  is  introduced  in 
considering  the  moments  acting  on  a  section  at 


r  ~ 


s 


0  <£■  s 


</)  1 

r*  -  r 


()■» 


The  origin  of  this  system  is  at  the  point  where  the  elastic  axis  intersects 
the  section.  Let  Q  ,  5  y(r>  ,  and  Q  yr)  be  the  components  of  the  total 
moment  acting  on  the  section  referred  to  these  axes.  {  5 x(r} and  3 y(r> 
are  taken  positive  in  the  positive  x(r) and  ^fr)- directions  but  5y,r,  is 
taken  positive  in  the  negative  y(r>- direction.)  It  then  follows  that 

«„r,  5  «*  '  <e"r  ^  pi  -  (eT*  >  P, 


Q 


yin 


i/i  - 


■t)p. 
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* sP* 


where  it  is  understood  that  4  ■  <?,  .  Py  .  etc. ,  are  force  and  moment  com¬ 
ponents  applied  at  station  *  r'f*  .while  rt  and  are  elastic  axis 
defections  at  station  r  .  When  second-order  terms  in  the  oscillatory 
variables  a^e  dropped,  these  equations  reduce  to  the  following; 


tin 


*■  ur. 
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*  et'3x  i  *  3  / 

V>  *  *yre'y  %  *  h  *  • 
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Here  the  steady  deflections  of  the  elastic  axis  have  been  included  in 
the  definitions  of  el ^  and  e^J  .  Since  the  blade  section  is  not  normal  to 
the  t  -axis  when  there  is  bending  deformation,  and  the  blade  chord  is  at 
an  angle  relative  to  the  j/(r;-axis,  the  above  moments  must  be 

resolved  accordingly.  'The  true  local  flatwise  bending,  chordwise 
bending  and  torsionaj.  moments  (which  are  to  be  equated  to  the  elastic 
restoring  moments  Mf,  M t  and  £  )  are  given  by 

=  ($$n)  +  "°5  </*+  -  (Q}tr>  ~  QX(r)^  ^ 


=  ®  ^  sm/^  f)  +  (Q}(rf  <f>) 


Q  =  Q  -  Q  i <■  Q  ,  ,  ur' 

*  *x(r)  *f<r)  ^ 

respectively.  After  dropping  small  terms  of  second  order,  substituting 
for  Q.tir)  .  dy(r;tand  Q  f(rj  ,  and  using  the  conversion  formulas  in 
Equation  (40),  the  above  equations  become 

ZZ  /S"  a  .  e  .  _  n  * 1  f'ry  \  1 


», 1  «. *•  *  K  K  *  *  pi  -  ?« ■ *,'  ■ *  <?.  - s  i  > * 
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where 
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Pk  -  Py  sm/S  ***■  Py  cosft 


Qa~  Qy  <-osj5t>"  Qy  sm / 


(*\  E  r„A(/> 


<?fc=  Qy  smfi*  t  Qy  cos  ft 
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p  -  P.  cos  3  -  Pu  5  m  £ ' 

a.  y  '  if  ' 


A  A  It)  £  „  </■> 

^  =  Pr  S'n/t*  rpy  cos/* 


^  =  (<3y  +  ef  Pj  « s/S<f)-($f  +  e^})  4  0 

Q  =  (<2y  +■  e'y  P,  )sin/(J,JT (q^+  Px)cosP 


0  s  Q  ,*•  eP  Pu  -  ej  P*.  . 

*  *  !  i)  J  f 


It  will  be  recalled  that 

2  w-  Ay  ,  urt  *  ur  +  Ly  , 

Ay  2  -j\y<P'cLX'  „  and 

Once  the  shape  of  *  vs.  x,  is  estimated,  say  ^if  ,  the  following  two 
quantities  can  be  computed  at  any  desired  r1  : 

£)s'»' ?dr l  e> ^  **m 

Er  er-*p 

Then,  at  the  desired  a:,  , 


Au  =  ^ 


A,  =  £,  * 


rt  -  r  -  Ey  <f> 


I 


\£r  * 


&1t  -  ur  cos /S  -  irt  smfi  =  &  t  *■  E z  4> 


'  where 


Szt  -  urf  smfi  +  vt  cosfi  =  &z~  £,<£ 
£  i  =  By  COS ft  -  E  y  s/n  /$ 

E  j  =  Ey  s/n +  E^  cos /S 


Since  the  steady  portions  of  the  moment  components  only 
cause  steady  displacements,  they  may  be  dropped  when  only  the  oscilla¬ 
tory  displacements  are  being  considered.  Thus,  the  following 
expressions  for  the  oscillatory  moment  components  are  obtained  from 
Equations  (47),  (48),  and  (49)  : 


»,  -  v  4  s,  - K  *(Sa-£, 


Q  -  Q  e'f)  P 

/«)  A  ^  *V  /  «v  / 

-  £  +  p.s<-paz,  +  s' 

<  A  b  •  *  b  "  <1  * 

+(p„s;-iK ). 

where 

Pa.-  P.-i 

pb -Pb-h  (P*m\-Pt 

m 

Q*.- 0«"  0* 

(§„■= «,«»/'’- «, 

0*  = 

(fl.  -  sin/l'v  -  cos/'’’) 

V 

(50) 


(51) 


(52) 
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*  ’ » 1 1  ^  b  ir  T  '  t  ti  r  i  •  -  tit'  ‘  i  >u  *  ,  \  4  -  • » *  ■  *  ^  j  ,>  i .  \  .  *  *  >  i »  1 1 .-  *  S »  « 

IV  h^fi  th^fr  a  !ol,ii  i.utubrr  mI  ^  *  nm  rtfif  fU  nut  ►  '%#-•*•  ouM^uffl  **1 

»"  ,  use  ui  the  turtnulas  in  Equation  (4b)  (wuh  wu,bein{«  the  ma«;.  ol 

the  ith  mass  bay  concentrated  at  r  -  r1')  yields 


o  u>  u> 

Pr  -  fl  2_  ^  r 
Z  <  =  / 


(V  -  1  for  the  tip  mass  bay. ) 


PA  ~  -(rCcos/3l/>)  (0£+/&r)/t,  *!C<) 

k  k 

-(SLX  s,n/*y)  '  £  mu,Cy  U)Aep  -  9A  rp  - /D  r,  )  Z,  V 


Ph  s  - (jft*  ^  (ee^^f)YL  in (,) r (° 


+  (fl*  cos/5(f))‘  iZ  -w  ulCjl)+  (ep  -0A  rp~/8B  r  )^T  *7f<; 

V 

■ 

3,  -  -  n‘  Z  («,  -v>  W,<V> 


+  (&£  +/3f  )H  nt,l)r  U>  cj‘ 


+  ■ri'ey>  *(‘>cJ'>t(ep~dA  rp  "A>  m  (t> 

*  ^<.-r  i-t 

-  (*,*/,)  £  mmru’ 

<L  -- -OlW*’)  ( £  -.'■V’cP -*•''£  me,rm) 

X  r  ‘‘I  '  7 

v'-’r'X’ -  e^’L  mu‘ru>) 


-(n‘„„/'’i£  rnV°CU‘  -  e.V,Z  w"'>rW, 

O  t;1  0  /  (sj 

-(n* cos/nit  rr,lorU)  Cf  -  e^L  mU)rU)) 

is  1  J  J  is} 
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.  i  t  j  '  •  *  i  I 


i  I  '  I  *. 


1 


The  ‘'ibration  analysis  is  carried  out  using  a  model  in  which 

the  blade  is  divided  into  spanwise  segments  or  bays,  each  of  which  is 
considered  to  be  a  uniform  beam  without  twist.  The  actual  twist  will 
be  taken  into  account  in  the  model  by  introducing  discontinuities  in 
twist  at  the  junction  between  bays  as  indicated  in  the  diagram  below. 


{A}“J=  [P]  {A}'""’ 


Figure  6.  ANALYTICAL  REPRESENTATION  OF  BLADE  TWIST. 
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I  i n  <  i  -  i  a  Ui».  i  uiut un  tti.il ici'-  |Af .which  tiff  di ;  phi v <'d 

*•«>  t  icuft-  t,  ,tr<  the  ilatuiM’  and  chordwise  iorce  components  (  Pa  .  Pb  ); 
th<  hendmg  moments  (  ,  <$b  );  the  torsional  moment  { <3X  ) ;  the  flatwise 

detleition  and  tlatuise  slope  (  5{  ,  5/  );  the  chordwise  deflection  and 
chordwise  slope  {  <52 ,  S2  );  the  torsional  deflection  <p  ;  and  the  flapping, 
lead-lag,  and  pitching  angles  (  6^ ,  ,  fip  ). 

Some  of  these  quantities  will  change  where  there  is  a  twist 
continuity  because  the  flatwise  and  chordwise  directions  are  defined 
ditferentiy  on  the  inboard  and  outboard  sides  of  the  junction.  The 
changes  in  {a}  across  the  twist  discontinuity  can  be  found  by  multiplying 
by  a  rotation  matrix  [r]  .  The  nonzero  elements  of  this  13  by  13  matrix 
are: 

R  {  ~  cos^il  b^-y3(olx)J  -  a,  ,  for  i  -  1  to  8  ; 


and 


R  -(  ~  I  .  for  -  9  to  1 3  . 


The  form  of  the[/?J  matrix  is  indicated  explicitly  on  Figure  6. 


ELASTIC  MATRIX 

The  1  3  by  1  3  elastic  matrixj_£j  relates  column  matrices  for 
the  inboard  and  outboard  ends  of  an  elastic  bay  of  span  £  according  to 

{*f  =  [£]  (a)  *r>  (53) 

where  {a?  is  the  column  matrix  indicated  in  Figure  6  .  Since  ,  Py  , 
<3t  ,  3f  ,  /£,  and  3 p  do  not  change  across  any  elastic  bay,  the  only 
nonzero  elements  in  the  1st,  5th,  9th,  11th,  12th  and  13th  rows  o£[z] 
are 

“•?  ^55  ~  ^99  ~  ^ 11, n~  ^ U.U  ^ ’ 3,13  * 
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Siina  c 


(>  f 


„  i"*  >  .  ,  ip, 

-  tPa 


and 


=  8*"  *  IP'“' 

i>  b  o 


there  are  two  nonzero  elements  in  both  the  2nd  and  6th  rows;  they  are 

£ 22  ~  £  t6  ~  1  and  ^ 21  ~  ^<5  =  ^  * 

The  rest  of  the  rows  of  elements  depend  on  the  elastic  properties  of  the 
bay  and  are  derived  in  the  following  paragraphs. 

The  differential  equations  for  the  deformations  of  an  elastic 
bay  are  obtained  by  equating  the  acting  moment  components  M t  ,  and  Q  , 
as  given  by  Equations  (50)  through  (52),  to  the  corresponding  restoring 
moment  components,  as  given  by  Equation  (18)  through  (20)  of  the  blade 
elastic  behavior  analysis.  In  the  latter  equations,  the  steady  terms  are 
dropped  (namely, in  m2  and  -r(i*-eteA)^in  Q  ).  With  5",  6"  ,  d '  and 
at'z  as  defined  in  Equation  (40),  T  =  PK, 

3  s  (£S2  -  et  EZz)/5' 


and 


A*  CrJ+  Px  4*  +  8ctJ'+(ZB1-2etEB2  +  e*  EI}  )(/8's* 


the  oscillatory  restoring  moment  components  may  be  written 
M,  =  £I1  (S"r  d'  p) 

Mz  =  EZ2  8"  ~(EIz  cl;  +  B)4' 


GL*Af‘-  3  8" 

rf 

--  c+'-  At  ** 


-A-  M 

El, 


51 


where 


es.)  (A 


/ 

* 


c  *  gJt  pt  (•  .(fa,  -  -j£ 

It  can  be  seen  that  when  the  expressions  for  Mt ,  Mz  ,  and  4  given  in 
Equations  (50)  and  (51)  are  substituted  into  the  above  equations,  three 
second-order  simultaneous  linear  differential  equations  for  61 ,  Sz  and 
<f>  are  obtained  in  which  some  of  the  coefficients  vary  with  s  .  * 

In  the  expressions  for  Mf ,  Mz  ,  and  (?  given  in  Equations  (50) 
through  (52),  the  terms  containing  6/  ,  <5*  ,  St  ,  S2  or  9  as  a  factor 
account  for  the  effects  of  the  variations  of  the  acting  moment  components 
along  the  bay  due  to  the  variations  of  the  deformations  along  the  bay. 
Thus,  by  keeping  the  bay  lengths  short,  the  variations  of  6f  ,  8,  ,  6t  ,  8t  , 
or  9  within  each  bay  may  be  assumed  to  be  linear  in  the  evaluation  of 


M,  ,  Mz  and  Q 

along 

the  bay. 

By  using 

8‘ 

“3 

,(C.«  ■) 

s“  ”) , 

6 .  S'"'’ 

-  8'  S  ; 

s" 

avg 

1 

=  l  ' 

’’-(S')0'’ 

}’  5'- 

and 

V 

-iC 

9 

4*  4*”- 

-  • 

the  expressions  for  M,  ,  M  2  and  <2  in  Equations  (50)  through  (52) 
become 

=  At  r  A3s  <-Aj  S* 

*  Z>,  ♦  Z)2s  -  D3s * 

<3  «  A,  +  XjS1 

Since  each  elastic  bay  is  considered  to  be  untwisted,  it  follows 
(within  each  bay)  that 

(&',  >  S',  or  <f>')  -  (d, .  S3>  or  <t>)  =  “  ^j(Si  »  »  or  P) 

and 

(Kor  Kh  zftrtef  or  5t)  =  -^r  (<?,  or  $2)  . 
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Consequently,  the  desired  differential  equations  are 

* 

v\  =  -  <iz  P  ♦-  -=  Y~  I'AjS1) 

S>  --  (<  r  /# 7)  537  *  D,s‘) 

-V*  *2^7  P-<  -J-  • 

Then,  by  subsequent  integrations  and  use  of  boundary  conditions  at 
r  *  r<ae>  or  3  =  0.  it  follows  that 

M  >«,•  *  f  v; 


*  lA j(A’  ~7**s  *t  *>s’) 


EIt  *  Z  ^ Is  "  i 

6,  =  5,*"'  -  (S2  fr>s  -(<*,’+  £§J-)(f#4s  -  <t>  K"s ) 

*  Jiij  &  *  T  S‘  5  *) 

*-  {*.* it*.- ■£** *  t 
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After  eliminating  (/V<*s  -  s  )  in  the  second  and  fourth  of  these  six 

equations  by  the  use  of  the  sixth,  the  first  five  yield  the  following  five 
equations  for  the  derivatives  of  the  deformations  averaged  from  s  =  0 
to  5  =  l  : 


When  the  expressions  for  At,  As  ,  Di  ,  2>}  ,  ,\  2  ,  and  A.,  are  substituted  into 
the  above  five  equations,  they  can  be  written  in  the  .form 
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where 


and  is  a  5  by  5  matrix  with  its  elements  as  listed  below: 


a. 


11  ~ 


/ 


a 


12  ~ 


1 

Z£l, 


ct 


1 3  S 


2£I, 


56 


57 


58 


By  substituting  the  expressions  for  A,  ,  Z>f  and  X1  in  the  expressions  for 
the  Fj’s  ,  the  latter  expressions  may  be  written 

{Pa, 

Q 

*: 

<5, 

p» 

*; 

S2 

** 

4 

where  Qc  j  is  a  5-x-iO  matrix  with  the  following  elements: 

c  _ _ IL 

"  2EI, 

c  -  —  ^ 

'2  "  EI1 
&ts  =  0 

C  -  -  ^  p 
~  £If  ** 


59 


61 


c 


59 


2_ 

C 


/  A  a  .  A  A 

T  P*~£1PZ^£1/i+£2Pb 


The  following  matrix  equation  is  obtained  by  combining  Equations  (54) 
and  (55): 


_  u 

-h  J 

-5/  £ 

Ay, 

-  a;  tv 

elv$ 

-  K  i 

-  <f'  l 

ray, 

V.  '  ^ 


=  [«]-  [C] 


(56) 
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Let[rf]  denote  the  5-x-lO  matrix  obtained  by  the  product 

M  =  H"  [c]  . 

Then  the  relations 


lead  to  the  following  expressions  for  the  nonzero  elements  in  the  3rd, 
4th,  7th,  8th,  and  10th  rows  of  the  elastic  matrix: 


£ 3h 

dih 

(h  =  1 ,  2,  and  5,  ...  10) 

£ll  * 

£*h  ' 

cC 

Zh 

(h  *  1,  2,3,  and  5.  £,  ...  10) 

£*;s 

/  *■  d 

'  2  4- 

t  7h 

*3H 

(h-  f,  2,  ...  6,  and  8,  9,  10) 

Er,Z 

U*,7 
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(h  -  1,  2,  7,  and  9,  10  ) 


£ah  '  a*h 


£88  =  1+**8 


£  10,  h  '  *5h 


£\0,lo"  *5,10  ' 


(h=  1,2,  ...<?) 


The  explicit  form  of  the  elastic  matrix [V] is  indicated  below: 


f>]  = 

1 

l  1 


d12 

***,> 

^ 15 

dft 

*n 

d'$ 

dl,,o 

d-ll 

^zz 

<2  2) 

It-  d1¥ 

^ Z5 

^Zi 

*21 

d2,,c 

i 

l  l 


d„ 

di) 

^3* 

dss 

dit 

Ud„ 

dJS 

d31 

d),tO 

*» 

*» 

dft 

U¥y 

1  *  d*e 

d¥f 

1 

d*J0 

dS2 

d>, 

dsi. 

dss 

d»t 

^57 

dSS 

dsi 

1 

i 
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MASS  MATRIX 


The  13  by  13  mass  matrix, [X) ,  relates  the  column  matrices 
tor  positions  just  inboard  and  just  outboard  of  an  assumed  concentrated 
mass,  such  that 


The  column  matrix  {a}  is  the  same  one  discussed  previously  which  has 
the  matrix  elements  shown  on  Figure  6  . 

Expressions  tor  the  elements  in  the  [rj  matrix  can  be  derived 
from  the  expressions  for  pA ,  pb  ,  ^ ,  and  <fx  which  are  given  by  Equations 

(41)  through  (45)  in  the  section  on  blade  inertia  loading.  When  Equa¬ 
tions  (41)  through  (45)  are  used  for  this  purpose,  the  mass  fn  appearing  in 
them  is  considered  to  be  the  total  mass  of  a  segment  of  the  blade  rather 
than  the  mass  per  unit  length.  These  equations  are  for  the  combined 
oscillatory  and  nonoscillatory  loadings;  but  recognizing  that^s  px  -  p*  , 
Pb  -  Pb~  >  etc.  .  ^  is  easy  to  separate  the  oscillatory  parts  of  the 
equations.  Steady  components  of  the  inertia  loadings  which  remain  in 
the  coefficients  of  the  oscillatory  variables  can  be  expressed  in  terms  of 
basic  mass  parameters  by  using  Equation  (46).  , 

Also, for  sinusoidal  oscillations  at  frequency  <*J  ,  the  second 
time  derivatives  of  the  oscillatory  variables  can  be  replaced  as  follows, 

Sit  = 

S2t  =  -  (s2-  E1  P)  (see  page  47) 

-  A}  or  '  ^ ($  > /* f  ’  Xd  ,  or  /p) 

When  Equations  (41)  through  (45)  have  been  modified  as  dis¬ 
cussed  above,  the  elements  of  the  [f]  matrix  are 


C  r 

'  tl 


13) 


r)(l  *  rrt  >w  1  *■  il  1  S/n  7 ,3  ' 

-  ~  ~  7r/fl  COS  5  S'O  J 

1 8  '  1 
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Fl,io  =  *  'Ic  («>1  +  n.XS'nz/)+  £2  F„-Ef  E18 


Fi.u  =  y«  ^\(r-rFt  6£  Cp  -/£  eDp)  cos  -  (eA  +/s)(Cp  smfi  +e2  +  qc 


-  771  Cl?  r  cosfi 


FVz  =  *  “>*  {  (r  -  rv)  s,n/1- epp  (&e  cosf-  &A  sm/S)-  eE  (ez  -  %  )  J 


+  77f fl* r p  s/nfl 


f.13-  W 


(r  -rp  )(0e  smjS +  &A  cosjS)  +-e2  +  rjc 


-  vrfL2  -|e p  cos/S-  (6£+/fr)(r- rp  )  smfi 


*  (et  rjc  )  cos  Zfi  -  ei  sm  Zfi 


Fz¥_  -  -  Trj  fi  r 

F„os-**1c  rtr  +  E2Fz+ 

F5¥.  =  F1S  -  WA*S"7/  COS/S 
F5s  -  77]  (u)  1  r  fL  COS  X jS  ’> 
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F98  =  -  TV  u)*e1  -  vrfl  *■  r\c  )  sin  2ft  +  ef  cos  2ft 
+  (e£+ftF)r  cos  ft  *  (ep  -rpeA-  r*  ft ^  )  smftj 

FV0  =  mU)1  r(ce2^~  "rAl  '  0»  3%  “  ^  K)  1c  COS/ 

+  (i*~  4* )  cos  2/5  i-  qc  (e?  cos  2^  -  ef  s/n  2ftij- 


+  Fz  Ffv.  "  £i  Ffe 


FV,  =  < 


(r'rF+e£cp-fii>ej>p)  ca 
-&**&)(«*  CpC)) 

-  ?n  fl  r  Cy 

F%n"  \tr-ro  +  e*f>eA)cy 

'  8£(k  +■  ev  Cy) 


+  cy 


Fv>  - 


> 

tvco1  |(r -rp)(9A  Cy  +  0E  Cy)  +  K 

Kz~  K1  *  CpC y  - (6£  +  ft £ )(r  -Pp)  Cy 


tv  n 
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The  reason  for  Eu  -  /  being  the  only  nonzero  element  in  the  3rd,  4th, 
7th,  8th,  and  10th  to  13th  rows  is  that  8{ ,  8f  ,  8^  ,  Sa  ,  <p  ,  ,  and 

/ 3p  do  not  change  across  a  concentrated  mass. 

The  explicit  form  of  the  [V]  matrix  is  indicated  below: 


[>]■ 
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METHOD  OF  SOLUTION 


When  the  blade  bays  are  designated  as  in  the  above  sketch, 
the  complete  matrix  equation  of  the  blade  is 


=  (  Jf  [£]“ V  [Rf  [ET°  \fT)  {Aj ^  .  (58) 


Inasmuch  as 

P  (r'W  =  dU,p)  -  P =  0 

( 

in  the  (Aj  column  matrix,  only  the  3rd,  4th,  7th,  8th, 

10th,  11th,  12th,  and  13th  columns  of  matrix  [F]u,are  carried.  Thus 
[Ej^is  multiplied  into  the  13-x-8  matrix  that  has  replaced  [f]'.  This 
gives  another  13-X-8  matrix  into  which  [/?J°is  multiplied.  The  process 
is  continued  until  the  last  inboard  matrix  is  reached.  As  a  result,  13 
equations  are  obtained  relating  the  variables  at  the  blade  root  and  the 
blade  tip.  The  last  three  equations  merely  state  that /3F  ,  J59 ,  and  jSp 
do  not  change  from  the  tip  to  the  root  and  need  not  be  retained.  The 
remaining  equations  can  be  written  in  the  following  matrix  form  where 
the  [«]  matrix  is  identical  to  the  first  10  rows  of  the  product  matrix  in 
Equation  (58). 
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In  addition  to  the  five  boundary  conditions  which  apply  in  all  cases, 
namely, 


(reef)  -  (.root) 
=  °2 


there  are  three  other  conditions  obtainable  by  considering  the  possible 
constraints  at  the  flapping,  lead-lag  and -pitching  hinges  of  the  blade. 
Take  the  pitching  hinge  as  an  example.  If  Ap  is  the  spring  constant 
(in  inch-pounds  per  radian)  of  the  pitch  control  mechanism,  then 


^  p  ftp  =  ^ 


+■ 


s,n  dA 
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or,  using  =  Q€  cos  *  Qi  sm/3  and  cos A~Qa  s"> A > 


.  f  Jr  oaf)  ^  Aroet)  >  (root) 

TCp/Sp  ~  (cos/s  s/n  QA  *■  smp  s/r\  0£  cos  QA) 


/  Jroot)  Jrott)  \  (root) 

+■  (sm JS  s/ndA  -  cos/S  sm&£  cosOA) 


+  ( cos  0£  cos  Qa  )  Q 


(root) 

% 


where 


„  (root) 


=  o( 


21 


(p;i 


(t,p) 


*  “zsAp 


C*‘} 


=  oc 


if 


*■  . 


*  «A» 


Q 


(roof) 

x 


oc 


1 1 


.  + 


“is  Ap 


Similar  equations  can  be  written  for  the  constraints  at  the  flapping  and 
lead-lag  hinges.  Thus,  eight  simultaneous  linear  homogeneous  equa¬ 
tions  in  the  eight  unknowns,  ...  ,fl  ,  and  fip  ,  can  be  obtained  from 

matrix  Equation  (59).  The  determinant  of  these  equations  is  plotted 
versus  various  trial  values  of  a)  for  given  rotor  angular  velocity,  A, 
and  nonoscillatory  flapping,  lead-lag,  and  pitching  angles  ( ,  and^ ) . 
The  zeros  of  this  curve  are  the  solutions  to  the  vibration  problem  for 
the  given  XI ,  /SF  ,  /Sa  ,  and  p?  .  The  simplest  case  is  the  one  with  all 
three  hinges  rigidly  locked  so  that/r*.  =  /fi  =  /SF  =  0  and  the  determinant 
is  simply 


oc 


31 


oc 


32 


OC 


Ji 


CC 


if 


OC 


OC 


ft 


’ll 


k St 


’tot 
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’103 


This  is  the  case  considered  in  Reference  2  . 


73 


EXPERIMENTAL  APPARATUS 


GENERAL  DESCRIPTION 


The  purpose  of  the  experimental  apparatus  is  to  apply  a 
controlled  vibratory  force  of  up  to  6  pounds  (12  pounds  peak  to  peak) 
to  each  blade  of  a  two-bladed  rotor  system  while  operating  in  a 
vacuum  at  rotational  speeds  up  to  1900  rpm.  The  equipment  which  has 
been  developed  has  the  following  capabilities.  The  vibratory  force  can 
be  applied  to  the  blades  so  that  either  the  symmetric  or  antisymmetric 
degrees  of  freedom  of  the  rotor  can  be  excited.  The  radius  at  which  this 
force  is  applied  is  adjustable  with  a  maximum  of  46  inches.  A  unique 
feature  of  this  equipment  is  that  the  angular  orientation  of  the  line  of 
action  of  the  vibratory  force  relative  to  the  plane  of  rotation  of  the  blades 
is  adjustable.  It  can  be  varied  from  lying  in  the  plane  of  rotation  to  lying 
normal  to  the  plane  of  rotation.  This  feature  enables  excitation  of  either 
the  chordwise  or  the  flapwise  degrees  of  freedom  and  also  allows  the 
input  generalized  force  in  a  specific  degree  of  freedom  to  be  maximized 
by  aligning  the  applied  force  with  the  motion  of  the  blade  at  the  point  of 
application.  The  vibratory  force  is  measured  at  the  attachment  point  of 
the  shaker  system  to  the  blades  with  a  specially  designed  force  coupling. 

The  test  apparatus  is  composed  of  two  major  assemblies  as 
installed  in  the  10 -foot-diameter  vacuum  tank  (Figure  8  ).  The  lower 
assembly  consists  of  the  test  rotor,  the  slipring  assembly,  the  hydraulic 
drive  motor,  and  the  base  supporting  structure.  The  upper  assembly 
consists  of  the  rotating  shaker  system,  a  50-pound  electromagnetic 
vibrator  and  the  overhead  supporting  structure.  These  two  assemblies 
are  joined  at  their  common  center  of  rotation  by  a  flexible  drive  coupling. 
This  coupling  not  only  transmits  the  driving  torque  to  the  upper  system 
but  also  fixes  the  azimuthal  alignment  of  the  rotating  shaker  system  rela¬ 
tive  to  the  test  rotor.  The  two  assemblies  are  held  erect  and  in  axial 
alignment  by  a  set  of  four  cables  and  turnbuckles.  The  hydraulic  motor, 
slipring  assembly,  and  electromagnetic  vibrator  are  in  sealed  enclosures 
to  allow  forced-air  cooling  while  the  system  is  operating  in  a  vacuum. 

The  driving  power  was  supplied  by  a  7-1/2-hp  electric  motor  and 
hydraulic  pump  mounted  outside  the  vacuum  tank. 

Four  energy-absorbing  safety  barriers  fabricated  of  wood  and 
metal  were  installed  so  as  to  surround  the  plane  of  rotation  of  the  blades 
and  the  shaker  system.  An  emergency  shutdown  device  was  provided 
which  monitored  the  vibration  level  of  the  supporting  structure  of  the 
lower  assembly.  If  the  vibrations  exceeded  a  preset  level,  the  device 
would  bring  the  rotating  system  to  a  stop  by  shutting  off  the  hydraulic 
pump  and  introducing  resistance  into  the  hydraulic  circuit. 

THE  ROTOR  HUB  AND  BLADES 


An  existing  rotor  hub  was  used  which  was  designed  and 
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fabricated  under  the  sponsorship  of  the  Air  Force  Aeronautical  Systems 
Division  for  use  in  wind-tunnel  tests  (Reference  1).  The  hub  can  be  seen 
in  Figures  8  and  9.  The  flapping  axis  for  each  blade  passes  through  the 
axis  of  rotation  and  normal  to  it.  The  lag  hinges  are  parallel  to  the  axis 
of  rotation  and  offset  5  inches  from  it.  Two  steel  flexure  straps  (3.  6 
inches  in  length)  connect  each  blade  to  its  lag  hinge  and  provide  a  pitch 
degree  of  freedom.  A  fixture  was  fabricated  to  interlock  the  two  sides 
of  the  hub  and  enable  it  to  be  converted  to  the  teetering  configuration 
when  desired. 


The  test  blades  are  uniform  beams.  They  were  fabricated  as 
solid  aluminum  bars  with  0.  5-x-3.  0-inch  rectangular  cross  section,  and 
when  installed  on  the  rotor  hub  they  have  a  tip  radius  of  44.  75  inches. 
They  were  mounted  so  that  the  effective  pitch  axis  of  the  rotor  hub  coin¬ 
cided  with  the  mid -chord  of  the  beams.  This  placed  the  blade  c.  g.  on 
the  pitch  axis  and  eliminated  the  pitch-flap  mass  coupling.  The  blade 
properties  for  the  uniform  section  outboard  of  the  retention  pins  are: 


W  eight 
Flapwise  El 
Chordwise  El 
Torsional  inertia 
Torsional  GJ 

Inboard  from  the  retention  pins  to 
the  pitch  flexure  straps  are: 

Flapwise  El 
Chordwise  El 
Torsional  GJ 

The  flapping  inertia  of  rotor  hub  p 
considered  infinitely  stiff  relative 


0.  1545  lbs. /in. 

0.  3125  x  106  lb.  -in.2 
0.  1125  x  108  lb.  -in.2 
0.  1155  lb.  -in.2 
0.436  x  106  lb.  -in.2 

s  lag  hinge, the  elastic  properties  of 

0.  1795  x  107  lb.  -in.2 
0.  529  x  106  lb.  -in.2 
0.  135  x  105  lb.  -in.2 

? 

blade  is  159  lb.  -in. ,  and  it  is 
the  blades. 


THE  LOWER  DRIVE  ASSEMBLY 

The  structure  of  the  lower  assembly  consists  of  two  large 
welded  box-like  enclosures  bolted  together  with  their  connections  sealed 
using  "0"  rings  and  rubber  gaskets.  The  base  section  encloses  the 
hydraulic  motor  and  has  removable  access  side  panels  with  bulkhead 
fittings  for  the  hydraulic  lines,  the  air  lines,  and  the  instrumentation 
leads.  The  top  section  encloses  the  slipring  assembly.  A  flexible 
coupling  connects  the  hydraulic  motor  to  the  rotor  shaft.  A  variable 
reluctance  pickup  was  installed  in  the  lower  enclosure  section  to  detect 
shaft  revolutions.  Inside  the  top  section  of  the  enclosure,  a  44-conductor 
slipring  was  mounted  on  the  shaft.  Above  the  slipring  assembly,  the 
shaft  was  hollow  with  outlets  on  either  side  of  the  top  bearing  and  seals. 
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SINGLE 


gure  9.  SHAKER  SYSTEM  AMD  ROTOR  BLADE  HUBS. 


Bulkhead  plugs  installed  on  these  outlets  allowed  the  instrumentation 
leads  on  the  rotating  side  o£  the  sliprings  to  be  passed  out  of  the  sealed 
enclosure  through  the  hollow  shaft  to  a  circular  terminal  board  below  the 
test  rotor  hub.  The  upper  end  of  the  rotor  shaft  is  splined  to  accept  the 
test  rotor  hub.  This  splined  section  of  the  shaft  is  also  threaded  to 
permit  vertical  adjustment  of  the  hub  position  on  the  shaft  by  means  of 
threaded  collars  on  either  side  of  the  hub. 

THE  ROTATING  SHAKER  SYSTEM  AND  UPPER  ASSEMBLY 

The  shaker  system  and  upper  assembly  are  suspended  over-  - 
head  on  the  center  line  of  rotation  of  the  rotor  system  by  a  truss 
structure.  The  bearing  assembly  for  the  rotating  shaker  system  and 
vibrator  mounting  are  attached  to  the  truss  by  four  threaded  rods  which 
allow  adjustment  of  the  vertical  position  of  the  shaker  system  and  angular 
adjustment  of  its  axis  of  rotation.  The  vibrator  is  mounted  above  the 
bearing  assembly  and  is  sealed  in  an  enclosure.  The  cooling  air  and 
electrical  lines  for  the  vibrator  pass  through  the  enclosure  via  bulkhead 
plugs  and  fittings. 

The  shaker  system  consists  of  a  rotating  hub  and  two  diamet¬ 
rically  opposed  torque  tubes  which  oscillate  about  their  longitudinal  axes 
(i.  e. ,  pitching  oscillation).  The  pitch  oscillation  of  the  torque  tubes  is 
driven  by  the  vibrator  through  a  series  of  push  rods,  rocker  arms  and  an 
inboard  pitch  arm  on  each  torque  tube.  Each  torque  tube  is  connected  to 
a  blade  through  an  outboard  pitch  arm,  force  rod  (push  rod),  and  force 
coupling.  Backlash  and  friction  are  minimized  by  the  use  of  flexures 
throughout  the  entire  oscillating  shaker  system  except  for  a  bearing  at 
the  connection  between  the  rotating  and  nonrotating  sides  of  the  system. 

The  advantage  of  this  shaking  system  is  the  capability  of 
changing  the  angular  orientation  of  the  shaking  force  relative  to  the  plane 
of  rotation.  The  line  of  action  of  the  shaking  force  can  be  varied  from 
no-rmal  to  the  plane  of  rotation  (as  when  forcing  the  blade  flapwise  modes) 
to  parallel  to  the  plane  of  rotation  (as  when  forcing  the  chordwise  modes). 
This  is  accomplished  by  changing  the  angular  orientation  (in  pitch)  of  the 
outboard  pitch  arms  and  force  rods  relative  to  the  torque  tube.  The 
radial  location  of  the  shaking  force  can  also  be  changed  by  changing  the 
radial  location  of  these  outboard  pitch  arms. 

The  installed  torque  tube  hub  can  be  seen  in  Figure  9.  This 
hub  rotates  in  and  is  carried  by  the  bearing  assembly  below  the  vibrator. 
A  push  rod  from  the  vibrator  (which  is  in  the  nonrotating  system)  extends 
down  through  the  hollow  shaft  of  the  hub  to  a  rotating  bearing  flexure 
assembly.  The  bearing  of  this  small  assembly  is  the  connection  between 
the  rotating  and  nonrotating  sides  of  the  shaker  system.  Short  links 
connect  this  bearing  assembly  to  each  of  the  rocker  arms  (Figure  10) 
which  are  mounted  on  the  top  surface  of  the  hub.  For  symmetric  shaking 
of  the  blades,  a  single  rocker  arm  is  used  for  each  torque  tube.  Anti¬ 
symmetric  shaking  is  obtained  by  replacing  one  of  the  single  rocker  arms 
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Figure  10.  ASSEMBLY  OF  SINGLE  ROCKER  ARMS  WITH  BEARING  FLEXURE. 


Figure  11.  DOUBLE  ROCKER  ARM  ASSEMBLY. 


with  a  double  rocker  arm  assembly  (Figure  11);  this  reverses  the  motion 
of  one  torque  tube.  A  push  rod  and  inner  pitch  arm  connect  each  torque 
tube  to  the  rocker  arms,  as  shown  in  Figure  12. 

Each  of  the  hollow  torque  tubes  is  cantilevered  from  the  hub 
by  two  pairs  of  pitch  bearings  which  react  the  transverse  loads  on  the 
torque  tubes.  The  two  torque  tubes  are  interconnected  by  a  tension 
strap  inside  the  hub.  The  axial  (centrifugal)  loads  of  the  torque  tubes 
are  thus  reacted  against  each  other  through  this  tension  strap. 

The  outboard  pitch  arm  and  force  rod  are  shown  in  Figure  13. 
The  pitch  arm  is  clamped  to  the  torque  tube  to  facilitate  both  radial  and 
angular  adjustments.  The  force  rod  is  connected  to  the  outboard  pitch 
arm  by  a  flat  flexure  designed  to  carry  the  centrifugal  shear  load  of  the 
force  rod.  The  force  coupling  which  transmits  and  measures  the  shaking 
force  is  assembled  at  the  lower  end  of  the  force  rod. 

This  force  coupling  accommodates  the  relative  angular  motions 
of  the  blade  with  respect  to  the  force  rod  without  significantly  influencing 
the  oscillating  shaking  force  it  transmits  and  measures.  The  complete 
assembly  and  an  exploded  view  of  its  components  are  presented  in  Figure 
14.  The  shaking  force  is  transmitted  through  a  pair  of  force  wires  pre- 
loaded  by  a  pair  of  small  beams  strain  gaged  to  measure'the  oscillating 
force.  The  centrifugal  load  of  the  coupling  is  reacted  by  the  blade 
through  a  tension  link.  A  small  aluminum  fixture  clamped  to  the  blade 
provides  an  adjustable  attachment  point  for  the  force  wires  and  tension 
link.  At  the  attachment  point  this  coupling  adds  0.  45  oz.  of  mass  to  the 
blade,  and  at  this  radial  location  it  adds  0.  82  oz.  of  tension  per  g  of 
centrifugal  acceleration.  These  small  additional  loads  are  included  in 
the  analysis,  and  their  effects  on  the  predicted  results  are  found  to  be 
small  for  the  test  blades  of  this  investigation. 
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•)  INSTALLED  VIEW 


b)  EXPLODED  VIEW 

Figure  14.  FORCE  COUPLING. 


82 


INSTRUMENTATION 


The  test  blade  was  instrumented  to  determine  the  spanwise 
distribution  of  moments  in  the  flapwise,  chordwise,  and  torsional 
bending  modes  by  measuring  the  moments  at  five  radial  stations.  These 
were,  nominally, r/R  -  0.  26,  0.  40,  0.  55,  0.  70,  and  0.  85.  As  a  check  on 
the  relative  phase  and  amplitude  of  the  response  of  the  two  blades,  the 
moments  were  measured  at  r/K=  0.  26  on  the  second  blade.  The  moments 
were  measured  using  C.  E.  C.  carrier-system  "D"  and  strain-gage 
bridges  with  four  active  arms. 

The  flap  and  lag  angles  of  the  test  blade  were  measured  with 
calibrated  strain-gage  beams  bridging  their  respective  hinges;  a  strain- 
gage  bridge  installed  on  the  pitch  flexure  straps  was  calibrated  to 
measure  the  blade  pitch  angle.  The  shaking  force  transmitted  to  each 
blade  was  measured  with  a  calibrated  four-active-arm,  strain-gage 
bridge  installed  on  the  force  beams  of  each  force  coupling.  The  rotor 
shaft  speed  was  determined  from  a  once-per-revolution  pulse  generated 
by  a  variable  reluctance  pickup  and  steel  pin  protruding  from  the  shaft. 
The  oscillating  input  voltage  to  the  electromagnetic  vibrator  was  recorded 
as  an  indication  of  the  total  force  input  to  the  shaker  system  and  the 
shaking  frequency. 

The  available  slipring  assembly  limited  the  rotating  instru¬ 
mentation  to  20  channels.  Twenty  C.  E.  C.  system  "D"  carrier 
amplifiers  and  two  power  supplies  were  used.  Four  power  rings  on  the 
slipring  assembly  allowed  the  instrumentation  to  be  commonly  powered  in 
two  groups  by  the  two  power  supplies.  The  data  were  recorded  on  two 
6-inch  (C.  E.  C.  tvpe  5-114)  oscillographs  which  were  slaved  together. 

An  internally  gene;  ad  correlation  trace  was  recorded  on  each,  and  an 
accurate  timing  rel-jrence  was  supplied  by  recording  a  60-Hz  signal  on 
each  oscillograph. 

In  addition  to  the  data  instrumentation,  nine  channels  of  strain- 
gage  instrumentation  were  required  for  monitoring  the  steady  and  oscil¬ 
lating  stresses  in  the  rotating  shaker  system  during  these  initial  tests. 
The  requirement  to  monitor  the  nine  channels  of  stresses  in  the  shaker 
system  left  only  eleven  channels  for  data.  The  eleven  data  signals 
recorded  simultaneously  were  the  three  blade  root  angles,  the  two 
shaking  forces,  five  of  the  test  blade  moments, and  one  moment  from  the 
second  blade.  The  six  moments  chosen  to  be  recorded  simultaneously 
were  either  the  flapwise,  chordwise,  or  torsional  moments,  depending  on 
the  mode  being  excited.  The  switching  of  signals  to  be  recorded  was 
facilitated  by  miniature  plug-connectors  installed  in  the  rotating  system. 
The  wiring  harnesses,  terminal  strips,  jumpers,  plug-connectors,  etc., 
in  the  rotating  system  were  arranged  to  provide  a  flexible  data  selection- 
capability,  e.  g.  ,  the  nine  channels  which  were  used  to  monitor  stresses 
during  these  tests  could  be  transferred  to  record  additional  blade-moment 
signals  by  changing  plug  connections. 
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The  50 -pound  electromagnetic  vibrator  was  driven  by  a 
variable  oscillator  and  power  supply.  A  frequency  counter  monitoring 
the  oscillator  f  -ovided  the  operator  with  an  indication  of  the  shaking 
frequency.  T1  j  operator  monitored  the  shaking  force  to  the  test  blade 
and  its  response  (the  moments  at  one  station)  on  the  x  and  y  axes  of  an 
oscilloscope.  The  rotational  speed  of  the  system  was  read  on  a 
frequency  meter. 
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EXPERIMENTAL  PROGRAM 


The  objective  of  the  experimental  phase  of  this  program  was 
to  develop  an  apparatus  and  a  testing  technique  for  determining  the  natural 
vibration  modes  and  frequencies  of  two-bladed  model  rotors  in  the 
absence  of  aerodynamic  loads.  The  initial  test  program  which  is 
described  herein  was  designed  to  determine  the  adequacy  of  the  apparatus 
and  testing  technique.  Thus  the  classical  configuration  of  uniform  blades, 
for  which  the  natural  vibration  modes  and  frequencies  can  be  computed 
accurately,  was  chosen  for  this  investigation. 

The  tests  were  conducted  in  a  vacuum  tank  with  the  air  density 
(and  thus  the  aerodynamic  loads)  reduced  to  about  1.  5  percent  of  that  at 
sea  level;  this  is  equivalent  to  an  altitude  of  approximately  100,  000  feet. 
The  blade  vibratory  modes  were  explored  at  rotor  speeds  of  300,  600, 

900  and  1200  revolutions  per  minute. 

It  is  noted  that  the  natural  vibration  modes  of  a  two-bladed 
teetering  rotor  include  all  the  modes  of  a  two-bladed  rotor  with  individual 
blade  flapping  and  zero  flap-hinge  offset.  This  assumes  that  the  blades 
and  hub  are  the  same  for  the  two  systems  and  that  the  rotor  shaft  imped¬ 
ance  as  seen  by  the  rotor  is  infinite.  The  rotor  hub  used  for  these  tests 
is  the  same  in  either  configuration.  It  is  converted  to  the  teetering  con¬ 
figuration  by  locking  together,  at  the  flap  hinge,  the  flapping  motion  of 
the  two  blades.  Therefore,  the  natural  vibration  modes  and  frequencies 
were  explored  in  the  teetering  configuration.  The  following  six  elastic 
natural  vibration  modes  were  investigated: 

1.  First  antisymmetric  pin-ended  flapwise  bending. 

2.  Second  antisymmetric  pin-ended  flapwise  bending. 

3.  First  symmetric  cantilever  flapwise  bending. 

4.  Second  symmetric  cantilever  flapwise  bending. 

5.  First  antisymmetric  pin-ended  chordwise  bending. 

6.  Second  cantilever  torsion. 

The  second  torsional  mode  was  investigated  instead  of  the  first  because 
the  pitch  flexure  straps  at  the  blade  root  were  torsionally  so  soft  rela¬ 
tive  to  the  blade  that  the  first  cantilever  torsional  mode  is  essentially  a 
rigid  body  pitch  response  with  virtually  no  torsional  deformation  of  the 
blade. 

Each  mode  was  investigated  by  bringing  the  rotor  up  to  the 
desired  rotational  speed  and  then  applying  the  appropriate  vibratory 
shaking  force  to  the  blades  (as  measured  by  the  force  couplings)  at  a 
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sequence  ot  closely  spaced  discrete  frequencies  over  a  frequency  range 
bracketing  the  natural  frequency.  An  oscillograph  recording  of  the 
structural  moments  and  shaking  force  was  obtained  at  each  shaking 
frequency.  The  shaking  force  and  the  structural  moments,  at  one  span- 
wise  station  on  the  blade,  were  monitored  on  the  x  and  y  axes  of  an 
oscilloscope.  This  facilitated  control  of  the  force  and  response 
amplitudes.  It  also  helped  to  nominally  locate  the  natural  frequencies 
by  observation  of  the  changing  Lissajous  pattern  due  to  the  phase  shift  of 
the  blade  response  relative  to  the  applied  force  as  the  shaking  frequency 
passed  through  the  natural  frequency.  The  natural  frequencies  are 
determined  by  the  mini  mums  in  the  shaking  force  per  unit  amplitude  of 
response  as  a  function  of  shaking  frequency. 
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THEORETICAL  COMPUTATIONS 


As  discussed  in  the  preceding  sections,  the  experimental 
investigation  of  the  natural  frequencies  and  mode  shapes  was  limited, 
in  this  particular  program,  to  a  rotor  system  with  uniform  blades. 

The  conditions  of  blade  uniformity  were,  therefore,  introduced  into 
the  general  theoretical  formulation  described  in  the  section  titled 
"Theoretical  Development,  "and  the  resulting  mass  and  elastic  matrices 
were  obtained  for  the  uniform  blade.  Two  existing  computer  programs 
using  the  matrix  method  of  solution  described  in  the  "Theoretical 
Development"  section  were  modified  to  use  these  mass  and  elastic 
matrices.  The  flapwise  and  chordwise  degrees  of  freedom  are  analyzed 
with  one  program  and  the  torsion  degree  of  freedom  with  the  other  pro¬ 
gram.  Both  programs  are  written  in  FORTRAN  IV  for  use  on  an  IBM 
360  computer. 

The  flapwise-chordwise  bending  program  consists  of  a  main 
program  and  8  subroutines.  It  has  a  storage  requirement  of  162,  000 
bytes.  The  blade  is  represented  by  a  series  of  concentrated  masses 
separated  by  massless  elastic  segments,  each  of  constant  stiffness. 

The  maximum  allowable  number  of  blade  segments  in  this  program  is 
100.  The  rotor  hub  is  represented  by  its  inertias.  The  program  com¬ 
putes  the  flapwise  or  the  chordwise  natural  frequencies  and  mode 
shapes;  it  also  has  the  capability  of  computing  the  coupled  flapwise- 
chordwise  modes  for  a  twisted  blade.  A  flow  chart  of  the  program  logic 
is  presented  in  Figure  15. 

The  segmentation  of  the  blade  representation  is  chosen 
to  represent  the  spanwise  blade  properties  adequately.  The  inputs 
are  the  spanwise  mass,  elastic,  and  twist  distributions  for  the  mode  to 
be  calculated;  the  hub  inertias;  root  boundary  conditions;  rotor  rota¬ 
tional  speed;  and  initial  trial  value  of  the  natural  frequency.  The 
printed  output  of  the  program  includes  all  the  inputs  in  addition  to  the 
computed  shear,  bending  moment,  slope,  and  deflection  at  the  spanwise 
position  of  each  concentrated  mass. 

The  computational  scheme  used  in  the  torsion  program  is 
based  on  the  Holtzer  method.  The  blade  is  represented  by  a  series  of 
concentrated  inertias  separated  by  inertialess  elastic  segments.  A 
flow  chart  of  the  program  is  presented  in  Figure  16. 
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BENDING  VIBRATION  COMPUTER  PROGRAM. 
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Figure  15.  (Cont'd) 


Figure  16.  FLOW  DIAGRAM  FOR  THE  TORSIONAL  BENDING  VIBRATION  PROGRAM 
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DISCUSSION  OF  TEST  RESULTS  AND 
COMPARISON  WITH  COMPUTED  RESULTS 


The  natural  frequencies,  mode  shapes,  and  corresponding 
shear  and  moment  distributions  were  computed  tor  each  of  the  six  blade 
modes  corresponding  to  those  ol  the  two-bladed  rotor  system  which  was 
investigated  experimentally.  Comparisons  o£  the  distributions  of  the 
measured  and  computed  spanwise  moments  and  the  natural  frequencies 
tor  each  of  these  modes  are  presented  in  Figures  17a  through  17£.  The 
comparisons  are  made  at  rotational  speeds  of  300  and  600  rpm  for  all 
modes  and,  in  addition,  at  900  rpm  for  the  first  symmetric  cantilever 
flapwise  bending  mode.  The  spanwise  moment  distributions  and  natural 
frequencies  computed  for  the  nonrotating  cases  are  also  presented  for 
reference,  although  they  were  not  determined  experimentally. 

The  computed  spanwise  moment  distribution  for  each  mode  is 
plotted  on  a  graph  as  a  continuous  curve,  and  the  measured  moments  are 
plotted  at  each  of  the  five  instrumented  stations.  The  amplitudes  of  the 
moment  traces  were  read  at  three  different  times  within  each  oscillo¬ 
graph  recording  {obtained  at  the  experimentally  determined  natural 
frequency),  and  the  resulting  three  sets  of  moment  data  are  shown  in 
the  figures.  For  identification,  a  different  plotting  symbol  is  used  to 
represent  each  of  these  three  sets  of  measured  moments.  The  spread 
of  the  measurements  at  each  spanwise  station  is  an  indication  of  the 
overall  accuracy  of  the  reduced  data. 

The  computed  moment  distributions  presented  are  for  the  mode 
shapes  normalized  to  a  unit  tip  deflection-  -  for  the  bending  modes,  one 
inch  of  tip  deflection  was  used;  for  the  torsion  mode,  one  radian  of  tip 
rotation  was  used.  The  measured  moments,  however,  were  obtained 
for  tip  response  amplitudes  considerably  less  than  these  unit  values. 
Therefore,  to  enable  comparisons  of  the  distributions  of  the  measured 
and  computed  moments,  the  measured  moment  distribution  in  each  mode 
was  scaled  so  that  the  maximum  measured  moment  was  made  equal  to 
the  computed  moment  at  the  same  radial  station. 

In  general,  the  measured  moment  distributions  are  considered 
to  be  in  satisfactory  agreement  with  the  computed  distributions,  but 
there  are  three  exceptions.  One  exception  is  the  antisymmetric  pin- 
ended  flapwise  bending  mode  at  600  rpm  (Figure  17a).  The  other  two 
exceptions  are  the  second  torsion  mode  at  300  rpm  and  at  600  rpm 
(Figure  17f). 

With  one  exception,  the  differences  between  the  measured  and 
computed  natural  frequencies  range  from  2  to  14  percent.  The  exception 
is  the  first  antisymmetric  p:n-ended  chordwise  bending  mode  for  which 
the  difference  between  the  computed  and  measured  natural  frequency  is 
about  35  percent. 
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Figure  17.  (Cont’d) 
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Figure  17. ( 
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MBS*  ANTISYMMETRIC  rid-iNDED  CH0R0WI S£  BENDING  MODE 
Figure  17.  (Cont’d)i 


Analyses  01  the  theoretical  and  experimental  results  have 
indicated  that  the  ditterences  can  be  attributed  to  both  experimental 
inaccuracies  and  limitations  in  the  assumed  mathematical  model. 

As  a  preeautionar\  measure,  the  shaking  force  during  these 
initial  tests  was  restricted  to  a  rather  low  level.  This  factor,  together 
cvith  the  high  stittness  ot  the  test  blades,  resulted  in  low  amplitude 
response  ot  the  blade  at  the  fundamental  shaking  trequeney.  Consequently, 
the  signal  lewis  were  \  erv  low,  in  general,  and  the  instrumentation  was 
operated  at  maximum  gain.  Internal  electronic  noise  was,  therefore, 
amplified  as  well.  Furthermore,  excitations  extrinsic  to  the  fundamental 
shaking  tori  e  were  apparentlv  being  introduced  to  the  system.  Sources 
ot  tlir  extrinsic  excitations  are  such  things  n»  the  rotating  shaker  system 
1 1 sell,  small  misalignments  ot  the  rotating  shaker  system  with  respect 
to  the  test  rotor  system,  the  hydraulic  drive  motor,  and  other  mechanical 
equipment  operating  in  the  building  (y  ia  the  vacuum  tank  structure).  These 
extraneous  excitations  and  the  ver^  loyy  structural  damping  of  the  solid 
test  blades  resulted  in  significant  responses  at  other  frequencies  in 
addition  to  the  fundamental  shaking  frequency  .  These  effects  superim- 
post-d  a  "muse "  on  the  signal  corresponding  to  the  response  to  the 
fundamental  shaking  torn*.  As  a  result  of  these  factors,  the  signal-to- 
noisi-  ley  els  on  the  records  yvere  quite  loyy  and  became  yvorse  yvith 
me  reusing  rotational  speed. 

Sample  oscillograph  records  are  presented  in  Figure  18. 

Both  ot  these  sample  records  yyere  obtained  for  the  first  symmetric 
cantilever  tlapwise  bending  mode  being  torced  very  close  to  its  natural 
t  requeue  y.  One  record  was  taken  at  a  rotational  speed  of  300  rpm  and 
tin-  othi  r  .it  '<00  rpm;  the  noise  level  and  its  increase  with  rotational 
speed  are  evident. 

Without  any  altering  ot  the  data  signals,  the  low  signal  >to- 
iioise  levels  resulted  in  readout  errors,  particularly  in  the  measurement 
ot  the  natural  frequent  n-s.  The  measured  natural  frequency  was 
determined  by  the  minimum  in  the  variation  with  fundamental  shaking 
trequeney  <d  the  -baking  tone  per  unit  amplitude  of  response  (i.  e.  ,  at 
the  peak  response).  Ihe  precise  definition  oi  this  minimum  was 
obscured  by  the  noise  level  in  the  data. 

Nevertheless,  while  there  are  limitations  on  the  experimental 
accuracy  lssoc  iated  with  problems  such  as  those  discussed  above ,  the 
probable  errors  m  the  measurements  are  not  sufficiently  large  to 
ii  i  omit  tor  the  dmerem  es  between  the  computed  and  measured  results 
m  general.  The  major  reason  tor  the  differences  that  have  been 
observed  m  tin-  investigation  is  that  the  analytical  representation  does 
not  correspond  to  the  actual  physical  system  that  was  tested.  It  was 
intended  that  the  model  configuration  and  conditions  designed  for  these 
initial  tests  should  bo  so  simple  as  to  be  representable  by  a  uniform 
*  km  m  having  negligible  intermodel  coupling.  This  objective  apparently 
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b.  900  rpm 


Figure  18.  (Cont'd) 
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was  not  achieved,  and  significant  effects  were  introduced  due  to  coupling 
and  the  simultaneous  excitation  of  several  degrees  of  freedom.  For 
example,  the  exceptionally  large  discrepancy  noted  previously  between 
the  measured  and  computed  natural  frequency  for  the  pin-ended  chord- 
wise  bending  mode  is  due  to  the  fact  that  the  measured  resonance  was 
actually  that  of  a  coupled-mode  response  whereas  the  analysis  assumes 
no  coupling.  Indications  are  that  the  measured  response  in  this  case 
was  a  coupled  mode  involving  not  only  pin-ended  chordwise  bending  but 
significant  amounts  of  torsion  and  motion  of  the  rotor  shaft  and  drive 
system  as  well.  In  addition  to  this  usual  type  of  linear  coupling,  there 
can  be  a  nonlinear  coupling,  such  as  a  change  (and/or  periodic  variation) 
of  the  generalized  mass  in  the  pitch  and  torsion  modes  proportional  to 
the  simultaneous  response  in  one  of  the  bending  modes.  This  intermodal 
coupling  due  to  multimode  forced  response  can  have  a  significant 
influence  on  the  measured  resonant  frequencies  of  the  system,  depending 
on  the  modes  and  their  amplitudes.  The  simple  analytical  representation 
which  was  assumed  does  not  account  for  the  response  of  the  rotor  on  its 
shaft  and  drive  system.  It  the  computations  for  the  pin-ended  modes,  the 
blade  moments  and  deflections  have  been  assumed  to  be  zero  at  the  pin; 
this  is  not  necessarily  true  for  the  test  rotor  system.  For  example,  the 
torsional  response  of  the  rotor  shaft  and  drive  system  will  allow  an 
oscillating  displacement  of  the  radially  offset  lead-lag  hinge.  Also,  the 
relatively  complex  tension-torsion  straps  used  in  the  hub  in  place  of 
pitch  bearings  and  a  pitch  spring  apparently  introduced  coupling  in  the 
system  which  is  not  accounted  for  in  the  analytical  representation. 
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CONCLUDING  REMARKS 


An  analysis  and  computational  technique  has  been  developed 
tor  predicting  the  natural  vibration  modes  and  frequencies  of  rotor 
blades  which  are  of  a  general  configuration  including  the  combined 
rigid  body  and  elastic  motions  of  the  blades.  This  was  programed 
for  the  special  case  of  the  uniform  blade  used  in  the  experimental 
phase  of  this  effort. 

An  experimental  apparatus  and  a  testing  technique  have  been 
developed  for  determining  the  natural  vibration  inodes  and  frequencies 
of  two-bladed  rotor  systems  in  all  degrees  of  freedom  while  operating 
in  the  absence  of  aerodynamic  loads.  Initial  tests  to  establish  the 
capabilities  of  this  experimental  apparatus  were  conducted  by  measuring 
the  natural  vibration  modes  and  frequencies  of  a  relatively  simple  two- 
bladed  rotor  system. 

The  experimental  phase  of  this  effort  has  demonstrated  the 
capability  of  the  experimental  apparatus  for  applying  controllable 
oscillating  shaking  forces  in  the  plane  of  rotation  or  normal  to  the  plane 
of  rotation  at  the  tip  of  each  blade  of  a  two-bladed  rotor  while  rotating 
in  a  near  vacuum.  This  program  has  also  demonstrated  the  feasibility 
of  the  experimental  technique  for  the  use  of  this  apparatus  to  determine 
the  natural  vibration  modes  and  frequencies  of  a  two-bladed  rotor 
system  while  rotating  in  a  near  vacuum. 

The  tests  performed  with  the  system  during  this  initial  pro¬ 
gram  have  indicated  that,  with  some  small  improvements,  sufficiently 
accurate  measurements  of  natural  vibration  modes  and  frequencies  can 
be  made.  Satisfactory  measurements  have  been  extracted  from  the 
experimental  data  even  though  the  signal-to-noise  levels  were  quite  low 
and  no  special  measures  were  taken  during  data  processing  to  minimize 
the  resultant  errors.  Some  elementary  mathematical  filtering  of  the 
data  might  have  improved  the  results  substantially. 

It  appears  that  the  extrinsic  excitations  constitute  the  primary 
limitation  on  the  adequacy  of  the  experimental  apparatus.  The  simul¬ 
taneous  excitation  of  the  rotor  at  frequencies  other  than  the  fundamental 
test  frequency  and  the  corresponding  response  can  have  a  significant 
influence  on  the  measured  resonant  frequencies.  The  coupling  effects 
(which  can  be  linear  or  nonlinear)  associated  with  such  simultaneous 
excitations  are  normally  neglected  but  could  be  important  depending  on 
the  particular  modes  and  their  amplitudes.  Some  effort  must  be  applied 
to  reducing  the  levels  of  these  undesirable  shaking  forces  or  to  isolating 
the  rotor  system  from  them. 

A  major  cause  of  the  differences  between  the  measured  and 
computed  natural  frequencies  appears  to  be  that  the  assumed 
mathematical  model  does  not  correspond  to  the  actual  system  tested. 
Some  of  the  complexities  in  the  test  model  were  introduced  as  a  result 
of  making  use  of  an  available  rotor  hub.  For  example,  the  pitch-tension 
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straps  used  in  this  hub  (in  place  of  a  pitch  bearing  and  pitch  spring)  are 
believed  to  have  produced  unexpectedly  large  coupling  effects.  However, 
the  model  rotor  system  used  in  these  tests  was  designed  on  the  basis  of 
standard  procedures  in  which  effects  such  as  those  described  above  and 
in  the  previous  section  have  normally  been  assumed  negligible.  Indica¬ 
tions  are  that  these  assumptions  are  not  valid  for  the  high-frequency 
modes.  Consequently,  although  these  initial  tests  were  intended  pri¬ 
marily  to  establish  the  capability  of  the  experimental  system,  they  have, 
in  fact,  already  demonstrated  the  inadequacy  of  structural  analysis 
methods  in  current  use. 


RECOMMENDATIONS 


!  In-  jjrrM'iit  ettort  has  developed  an  analysis  capable  of 
,ir«H!u  t.ng  the  d\namic  characteristics  of  rotor  systems  of  a  very 
general  «  mu .g>. rat ion  and  has  developed,  fabricated,  and  performed 
initial  lot*  with  the  experimental  apparatus  required  to  verify  the 
uuuUms,  While  the  teuMl>iht\  ot  using  the  developed  apparatus  and 
test  technique  hts  been  demonstrated,  there  is  a  problem  with 
extrinsic  excitation  ot  the  test  rotor  which  must  be  attenuated.  Also, 
because  ot  the  need  tor  \er:t\ing  the  analytical  procedure  which  has 
been  developed,  it  is  recommended  that  the  test  apparatus  and  technique 
<t»w  eloped  m  this  program  be  used  to  obtain  the  validating  data. 

Spew  itu  allv,  it  i»  recommended  that: 

1.  An  ettort  should  be  made  to  attenuate  the  spurious 
excitation  reaching  the  test  rotor  by  isolating  the 
experimental  apparatus  from  the  vacuum  tank  and 
trom  the  hydraulic  drive  motor. 

1.  For  anv  turther  tests,  relays  she-,  i  1  be  incorporated 
on  the  rotating  side  of  the  slipring  assembly  to  switch 
trom  stress-monitoring  signals  to  data  signals.  This 
would  permit  simultaneous  recording  of  flapwise, 
chordwise,  and  torsional  moments. 

}.  The  analytical  technique  developed  for  the  general 
rotor  configuration  should  be  programed  for  digital 
c  ompidation. 

-1.  Several  model  rotor  configurations  should  be  designed, 
lubricated,  and  tested  to  obtain  the  experimental  data 
necessarc  to  verity  the  general  analytical  method  that 
has  been  developed. 
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IS.  SNONSORINE  MILITARY  ACTIVITY 

U,  S.  Army  Aviation  Materiel 
Laboratories,  Fort  Eustis,  Virginia 

An  analysis  and  computational  procedure  has  been  developed  for  predicting 
the  natural  vibration  modes  and  frequencies  of  rotor  blade  and  hub  systems.  This 
rrethod  is  capable  of  treating  rotor  systems  which  are  of  a  general  configuration  and 
o'  including  configuration  variables  which  are  normally  neglected. 

An  experimental  apparatus  and  a  testing  technique  have  been  developed  for 
obtaining  the  experimental  data  necessary  to  verify  the  analytical  procedure.  This 
apparatus  and  the  technique  are  capable  of  determining  the  natural  vibration  modes 
and  frequencies  of  two-bladed  rotor  systems  of  a  very  general  configuration  while 
operating  in  the  absence  of  aerodynamic  loads.  An  initial  test  program  with  the 
apparatus  was  conducted  to  determine  the  adequacy  of  the  equipment.  The  tests  not 
only  demonstrated  the  capability  of  the  equipment  and  test  procedure  but  also 
revealed  a  problem  with  extrinsic  excitation  of  the  test  rotor  at  frequencies  other 
than  the  shaking  frequency.  Recommendations  are  made  for  modifications  to  the 
system  and  the  techniques  in  order  to  reduce  the  levels  of  these  extraneous  excita¬ 
tions  and  to  minimize  their  effects. 

The  computed  and  measured  spanwise  moment  distributions  and  the  corre¬ 
sponding  natural  frequencies  are  presented  for  six  of  the  natural  vibration  modes  of 
the  test  rotor  at  several  rotational  speeds.  The  theoretical  and  experimental  results 
are  compared  and  their  differences  are  discussed. 
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